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Abstract 


The  deformations  that  occur  around  shafts  and  tunnels 
depend,  in  part,  on  the  elastic  deformation  moduli.  These 
moduli  are  usually  stress-path  dependent  and  it  is  therefore 
important  to  follow  the  correct  stress  path  in  the 
laboratory  when  determining  deformation  moduli. 

The  stress-path  dependency  of  Athabasca  oi 1  sand  has 
been  investigated.  Three  series  of  drained  triaxial  tests 
were  conducted  on  'undisturbed'  samples  of  oil  sand  along 
different  stress  paths.  Two  of  the  series  were  conducted  in 
a  high  pressure  triaxial  cell,  and  the  other  series  was 
conducted  on  hollow  cylinder  samples  under  plane-strain 
conditions.  The  laboratory  results  have  shown  that  the 
deformation  of  oil  sand  is  path  dependent.  Of  even  greater 
significance  is  the  fact  that  sample  preparation  techniques 
have  a  dominant  role  in  the  sample  disturbance  that  occurs. 

A  numerical  model  was  developed  to  predict  the 
deformations  that  occur  in  the  hollow  cylinder  tests.  The 
predicted  deformations  are  good,  provided  that  the  Young's 
modulus  is  correctly  specified  and  that  the  model  accounts 
for  the  variation  of  this  modulus  with  confining  pressure. 
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CHAPTER  1 


INTRODUCTION 


1 . 1  Genera  1 

The  oil  sands  of  the  Athabasca  deposit  in  Alberta  are 
currently  recovered  by  open  pit  mining  techniques.  There 
are  now  two  such  mines  operating  near  Fort  McMurray  ; 
Syncrude  Canada  Ltd.  and  Suncor  Inc.  Two  further  operations 
are  scheduled  to  begin  construction  in  the  near  future  ; 
Alsands  near  Fort  MacKay  (open  pit  mining)  and  the  Cold  Lake 
Project  in  Eastern  Alberta  (in  situ  recovery). 

Open  pit  mining  is  only  feasible  for  deposits  that  do 
not  lie  deeper  than  approximately  70  m.  Much  of  the  rich 
oil  sands  in  Alberta  occur  in  the  middle  McMurray  Formation 
which  in  most  places  far  exceeds  the  practical  limits  for 
open  pit  mining.  As  the  world  price  of  petroleum  continues 
to  escalate,  it  is  now  becoming  economical  to  proceed  to  in 
situ  recovery  processes  which  involve  shafting  and 
tunnel  ling. 

Three  in  situ  recovery  processes  involving  tunnelling 
are  possible  (Devenney  and  Rai sbeck , 1 980 ) .  The  first 
involves  sinking  a  shaft  through  the  overlying  formation 
into  the  limestone  which  lies  beneath  the  oil  sands,  and 
then  tunnelling  into  the  limestone.  Injection  wells  would 
be  drilled  from  the  tunnel  into  the  oil  sand  above.  Another 
possibility  is  that  used  in  the  USSR  at  Yarega  (Devenney  and 
Raisbeck, 1980) .  The  shaft  is  sunk  to  a  depth  just  above  the 
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reservoir  and  a  tunnel  driven  horizontally.  Injection  wells 
are  drilled  down  into  the  reservoir.  The  third  arrangement 
is  similar  to  the  first,  the  tunnel  being  driven  in  the  oil 
sand  deposit.  Figure  1.1  shows  these  three  possible 
recovery  procedures . 

Irrespective  of  which  procedure  is  used,  it  is 
important  to  determine  how  the  oil  sand  will  behave  during 
excavation.  Shafting  and  tunnelling  in  either  the 
overburden  or  the  limestone  bedrock  pose  convention  problems 
only  as  the  behaviour  of  these  materials  is  well  understood. 
The  oil  sands  pose  special  problems  as  its  behaviour  upon 
excavation  is  not  well  understood  so  far.  The  problem  is 
further  complicated  by  the  lack  of  any  previous  experience 
of  deep  excavations  in  this  material. 

There  are  only  two  case  histories  available  to  be 
studied  and  neither  of  these  involved  deep  excavations.  The 
GCOS  test  shaft  (Hardy  and  Scott, 1978)  was  sunk  to  only  23.5 
metres  when  instability  of  the  walls  caused  severe  slabbing 
and  thus  the  abandonment  of  the  shaft.  The  Saline  Creek 
Tunnel  (Smith  et  a/., 1978)  showed  the  oil  sands  at  that 
location  to  behave  competently,  but  the  overburden  at  this 
site  was  only  27.5  meters. 

It  is  thus  imperative  that  more  geotechnical  knowledge 
be  obtained  and  used  with  all  the  available  theory  of  oil 
sand  behaviour  to  enable  a  better  understanding  of  how  this 
material  will  behave  when  excavated  to  depths  exceeding  300 


metres . 
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Figure  1.1  Possible  in  situ  recovery  methods  for  oil  sands 

(after  Devenney  and  Rai sbecK, 1980 ) 
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The  oil  sands  in  situ  are  a  very  dense,  uncemented 
sandstone.  Their  shear  strengths  are  much  higher  than  those 
of  normal  granular  deposits  due  to  the  interlocking  nature 
of  the  sand  grains  which  results  in  a  material  that  has  a 
very  high  dilatancy  ( Dusseau 1 t , 1 977 ) .  There  is  no  chemical 
cementing  present  in  the  oil  sands  and  this  results  in  a 
Mohr-Coulomb  peak  envelope  which  passes  through  the  origin 
(Dusseault  and  Morgenstern , 1 978 ) ,  and  hence  the  oil  sands 
may  be  classified  as  a  cohesionless  material.  At  higher 
stress  levels,  dilatancy  will  be  suppressed  and  this  results 
in  a  curvature  of  the  Mohr  envelope. 

A  tangent  drawn  to  the  envelope  (at  high  stress  levels) 
will  give  a  cohesion  intercept  which  has  been  defined  as  the 
fabric  cohesion  ( Dusseau 1 t , 1 977 ) . 

The  pores  of  the  oil  sands  are  filled  with  water, 
bitumen  and  gas,  which  is  present  either  as  free  gas  or 
dissolved  gas  (Dusseault  and  Morgenstern , 1 978 ) .  The  bitumen 
does  not  contribute  to  the  strength  of  the  material  at 
normal  temperatures,  but  has  a  marked  effect  on  the 
permeability  of  water  through  the  material.  The  oil  sands 
are  thus  low  permeability  materials  and  behave  in  an 
undrained  manner  when  subjected  to  a  fast  stress  or 
temperature  change  ( Dusseau 1 t , 1 979 ) . 

These  phenomena  of  oil  sands  are  obviously  factors 
which  would  contribute  to  the  behaviour  of  the  material  upon 
excavation.  When  the  material  is  unloaded  (or  the 
temperature  increased)  the  gas  comes  out  of  solution  and, 
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because  of  its  low  permeability,  there  is  an  expansion  of 
the  material.  This  expansion  results  in  a  disruption  of  the 
fabric,  and  hence  a  decrease  in  the  shear  strength  (Harris 
and  Sobkowicz, 1977 ) . 

The  exsolution  of  gas  continues  as  the  unloading 
progresses.  The  gas  bubbles  eventually  form  a  connected 
drainage  system  and  bleed  off  of  the  gases  occurs 
( Scot t , 1 978 ) .  At  this  stage,  the  material  begins  to  behave 
in  a  drained  manner.  One  of  the  major  problems  yet  to  be 
solved  is  the  question  of  how  long  a  time  interval  is 
required  before  drainage  occurs.  This  problem  is  of 
significance  in  the  design  of  shaft  and  tunnel  linings  and, 
until  solved,  the  behaviour  of  oil  sands  upon  excavation  can 
only  be  estimated  as  lying  between  its  two  modes  of 
behaviour;  undrained  and  drained. 


1.2  Aim  of  the  Thesis 

The  aim  of  this  thesis  is  to  examine  the  stress-strain 
behaviour  typical  for,  and  applicable  to,  ground  response 
around  tunnels  and  shafts.  The  presence  of  gas  will  have  a 
major  influence  on  the  undrained  behaviour  of  oil  sands  and 
it  is  thus  necessary  to  obtain  samples  for  laboratory 
testing  that  still  have  gas  in  the  material.  This  may  be 
accomplished  if  in  situ  freezing  techniques  are  used,  and 
the  problems  of  obtaining  undisturbed  samples  by  this 
technique  have  been  discussed  by  Dusseault  (1980).  The  gas 
has  no  effect  on  the  results  of  drained  tests  provided  the 
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samples  have  not  been  disturbed. 

The  research  undertaken  by  the  author  was  restricted  to 
examining  the  drained  behaviour  of  the  oil  sands  at  depths 
exceeding  300  metres.  The  tests  were  conducted  on 
relatively  undisturbed  samples  in  order  to  examine  the 
stress-deformation  behaviour  of  the  in  situ  material.  The 
data  obtained  from  such  tests  allows  for  an  evaluation  to  be 
made  as  to  whether  the  material  in  its  undisturbed  state  is 
competent  enough  to  withstand  the  stresses  imposed  upon  it 
during  shafting  or  tunnelling.  This  in  itself  is  only  one 
of  the  'bounds'  that  need  to  be  examined  to  fully  understand 
the  behaviour  of  oil  sands,  the  undrained  behaviour  being 
the  other  bound. 

The  laboratory  program  consisted  of  performing  a  series 
of  drained  hollow  cylinder  triaxial  tests  as  it  is  believed 
that  the  stress  path  followed  in  this  test  more  closely 
resembles  that  which  occurs  in  the  field.  The  influence  of 
stress  path  on  the  deformation  behaviour  was  examined  by 
conducting  two  other  series  of  triaxial  tests  on  solid 
samples  of  37.5  mm  diameter.  The  test  results  from  all 
tests  can  then  be  compared. 

The  analysis  of  the  hollow  cylinder  triaxial  data 
depends  upon  the  assumptions  made  regarding  the  stress 
distribution  within  the  cylinder.  A  numerical  model  was 
developed  to  interpret  the  test  results  and  thus  enable  some 
conclusions  to  be  drawn  concerning  the  validity  of  these 
assumpt i ons . 
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1.3  Outline  of  the  Thesis 

The  sample  collection  techniques  and  the  unique  method 
of  sample  preparation  for  both  the  hollow  cylinder  and 
triaxial  tests  are  discussed  in  Chapter  2.  Laboratory 
techniques  and  equipment  are  also  presented  in  this  chapter. 

Chapter  3  deals  with  the  theoretical  stress 
distribution  within  the  hollow  cylinder,  and  the  method  of 
calculating  strains  from  the  volume  changes  measured  during 
the  test.  In  Chapter  4,  the  development  and  theoretical 
aspects  of  the  numerical  model  are  presented.  A  sensitivity 
analysis  of  the  influence  of  the  various  parameters  on  the 
solution  is  given. 

The  laboratory  results  are  given  in  Chapter  5  along 
with  a  discussion  of  these  results.  The  hollow  cylinder 
triaxial  results  are  analysed  in  more  detail  in  Chapter  6 
which  compares  the  actual  test  results  with  those  predicted 
by  the  numerical  model.  A  detailed  discussion  of  these 
comparisons  is  given. 

The  conclusions  reached  from  this  research  are 


summarized  in  Chapter  7. 


CHAPTER  2 


LABORATORY  EQUIPMENT.  SAMPLES  AND  PROCEDURES 
2.1  Introduction 

The  stress-deformation  behaviour  of  oil  sands  is  being 
investigated  to  determine  the  behaviour  of  this  material 
when  deep  excavations  such  as  shafts  and  tunnels  are 
constructed.  It  was  decided  to  use  a  hollow  cylinder 
triaxial  device  (HCTD)  for  the  testing  program  as  this 
device  is  more  versatile  than  a  conventional  triaxial  cell 
and  allows  several  stress  paths  to  be  followed  in  loading  of 
the  sample.  A  stress  path  resembling  that  which  occurs  in 
the  field  can  also  be  modelled  in  this  apparatus  by 
maintaining  a  constant  cell  pressure  and  decreasing  the 
internal  pressure. 

The  HCTD  used  in  this  research  program  was  designed  and 
built  at  the  University  of  Alberta  ( Dusseau 1 t , 1 979 ) ,  and  has 
the  capability  of  sustaining  cell  pressures  of  up  to  20  MPa. 
These  high  cell  pressures  are  of  the  same  order  as  the  in 
situ  stresses  that  will  be  encountered  in  deep  excavations. 
Discussion  of  this  apparatus  and  the  problems  encountered 
will  be  dealt  with  in  a  later  section. 


2.2  Literature  review  of  previous  HCTD  usage 

Hollow  cylinder  tests  in  the  laboratory  have  been 
conducted  as  far  back  as  1957.  The  author  has  undertaken  a 
literature  review  to  determine  the  various  capabilities  of 
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the  different  HCTD' s  and  the  purposes  for  which  they  were 
employed.  This  review  has  been  limited  mainly  to  the 
testing  of  soils  since  the  uncemented  nature  of  the  oil  sand 
makes  this  material  similar  to  a  dense  sand,  rather  than  a 
rock . 

The  majority  of  the  researchers  have  used  the  HCTD  to 
investigate  the  shape  of  the  yield  surface  for  their 
particular  material.  The  HCTD  is  well  suited  for  this  type 
of  investigation  as  control  of  all  three  principal  stresses 
is  possible.  Thus  the  shape  of  the  failure  surface  in  3-D 
stress  space  between  the  two  conventional  points  (triaxial 
compression  and  extension)  may  be  obtained  by  following 
various  stress  paths  during  the  testing  program. 

Ki rkpatr i ck ( 1 957 )  used  a  medium  to  fine-grained  sand  and 
allowed  full  drainage  during  shear.  He  conducted  an 
extensive  testing  program  as  did  Wu  et  a/. (1963)  to 
determine  the  applicability  of  the  various  failure  criteria. 
Wu  et  al .  (1963)  conducted  their  experiments  on  both  a 
remoulded  plastic  clay  and  a  uniform  Ottawa  sand.  No 
drainage  was  allowed  during  shear,  but  pore-pressure 
measurements  were  recorded. 

Haythornthwai te( 1960 )  investigated  the  yield  and  flow 
criteria  of  a  remoulded  silt  and  his  apparatus  allowed  the 
application  of  both  a  vertical  stress  and  torque.  Frydman 
et  a/. (1973)  used  glass  microspheres  as  their  material  and 
ran  drained  tests  in  their  program. 

The  only  investigators  to  have  used  'undisturbed' 
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samples  in  an  experimental  program  were  SuKlje  and 
Dr novseK ( 1 965 ) .  They  used  a  Tertiary  clay  and  ran  a  series 
of  plane-strain  and  plane-stress  tests  to  determine  the 
tensile  deformabi 1 i ty  of  this  material.  Their  apparatus  did 
not  have  the  capability  of  allowing  drainage  nor  to  measure 
pore-pressure .  Details  of  their  sample  preparation  for  the 
undisturbed  clay  were  not  given.  Whitman  and  Luscher ( 1 962 ) 
ran  a  series  of  experiments  to  investigate  the 
soi 1 -structure  interaction  problem  of  a  buried  pipe.  The 
tests  were  conducted  on  dry  Ottawa  sand  under  plane-strain 
condi t ions . 

Broms  and  Jamal (1965)  used  the  HCTD  to  analyse  the 
results  of  conventional  triaxial  tests.  They  ran  undrained 
tests  on  a  rounded  beach  sand,  and  used  the  hollow  cylinder 
data  to  explain  why  the  angle  of  friction  measured  in 
triaxial  compression  and  extension  differed.  Their  results 
show  that  the  difference  is  due  to  a  non-uniform  stress 
distribution  in  the  sample.  In  another  paper  by  Broms  and 
Casbar i an ( 1 965 ) ,  the  hollow  cylinder  data  was  used  to 
examine  the  effects  of  rotation  of  the  principle  stress  axes 
on  shear  strength.  Procter ( 1 967 )  used  the  hollgw  cylinder 
to  investigate  the  stress-di 1 atancy  behaviour  of  a  dense 
sand . 

The  capabilities  of  these  various  devices  are  very 
similar.  All  authors  had  independent  control  of  the  cell 
pressure,  bore  pressure  and  axial  stress.  Volume  changes  in 
the  bore  could  also  be  measured,  while  the  volume  change  of 
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the  pore  fluid  was  only  measured  by  Ki rkpatr ick( 1 957 )  and 
Frydman  et  a/. (1973).  The  other  authors,  except  SuKlje  and 
Dr novsek ( 1 965 ) ,  had  the  capability  of  measuring  pore 
pressures .  Axial  strains  and  stresses  were  measured  by 
various  devices. 

The  cell  pressure  used  by  these  authors  all  varied,  but 
in  no  instance  did  it  exceed  400  KPa.  This  pressure  is  far 
less  than  the  proposed  pressures  to  be  used  in  this  testing 
program.  However,  the  use  of  these  lower  pressures  enabled 
the  researchers  to  use  a  modified  triaxial  cell  whereas  this 
research  requires  the  use  of  a  completely  new  piece  of 
equ i pmen  t . 

The  stress  distributions  within  the  hollow  cylinder 
specimen  were  calculated  in  various  ways.  SuKlje  and 
Drnovsek( 1 965 )  based  their  solution  on  the  theory  of 
elasticity.  They  realized  that  these  solutions  were  only 
valid  in  a  limited  stress  range.  Broms  and  Jamal (1965) 
conducted  triaxial  compression  tests  on  their  cylinders  and 
thus  had  a  uniform  stress  distribution  within  the  cylinder. 

The  research  carried  out  by  Frydman  et  al .  (1973), 
Whitman  and  Luscher ( 1 962 )  and  K i rkpatr i cK ( 1 957 )  was  only 
concerned  with  the  stresses  at  failure,  and  an  average 
tangential  stress  was  calculated  from  equilibrium 
requirements  of  a  half  section  of  the  cylinder,  the  average 
radial  stress  being  found  either  from  the  failure  criterion 
used  or  by  assuming  a  linear  distribution  of  the  radial 
stresses  within  the  cylinder.  Wu  et  al .  (1963)  determined 
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the  actual  stress  distribution  at  failure  from  a  limiting 
equilibrium  analysis. 

Procter ( 1 967 )  discussed  the  means  available  for 
calculating  the  stress  distribution.  He  decided  on  defining 
a  radial  stress  distribution  (linear)  as  he  did  not  think 
that  elaborate  elasto-plastic  analysis  was  valid  for 
particulate  media.  His  analysis  was  based  on  average 
stresses  since  the  calculated  strains  in  the  cylinder  were 
average  strains. 

Stress  distributions  can  also  be  found  in  text  books 
such  as  Jaeger  and  Cook(1976)  and  Timoshenko  and 
Goodi er ( 1 970 )  for  hollow  cylinders.  Jaeger  and  Cook(1976) 
give  solutions  based  on  elasticity  and  plasticity  while 
Timoshenko  and  Goodi er ( 1 970 )  only  deal  with  elasticity. 

The  literature  review  revealed  valuable  insight  into 
the  problems  associated  with  hollow  cylinder  testing. 
However,  no  information  was  found  that  would  assist  in  the 
sample  preparation  of  undisturbed  material,  and  the  low 
pressure  range  in  which  these  hollow  cylinder  devices  were 
operated  did  not  provide  much  insight  into  the  problems  that 
would  be  encountered  when  testing  at  high  pressure. 


2.3  The  University  of  Alberta  HCTD 

This  apparatus  is  described  in  detail  in  a  paper 
published  by  Dusseaul t ( 1 98 1 ) ,  a  copy  of  which  is  included  in 
Appendix  C.  The  basic  features  described  in  this  paper  have 
not  been  changed,  but  a  few  minor  modifications  were 
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necessary. 

The  top  cap  was  modified  such  that  the  porous  stone  can 
be  attached  to  it  and  be  flush  with  the  bottom  of  the  cap. 
This  eliminated  the  possibility  of  puncturing  the  external 
membrane  against  the  porous  stone  when  applying  high  cell 
pressures.  A  problem  still  existed  with  the  internal 
membrane  as  it  tended  to  pinch  between  the  stone  and  cap. 
This  problem  was  finally  solved  by  inserting  a  ring  of 
Neoprene  rubber  around  such  joints. 

It  was  desireable  to  have  the  capability  of  deairing 
the  system  after  it  was  assembled.  Two  additional  ports 
were  drilled  through  the  top  cap  such  that  they  opened  onto 
the  porous  stone.  Deairing  is  thus  possible  by  applying  a 
back  pressure  to  the  pore  pressure  line  and  opening  the  line 
at  the  top  of  the  cap. 

The  major  modifications  made  were  to  the  membranes. 
Several  leaks  in  the  internal  membrane  occurred  during 
testing  and  it  was  thus  necessary  to  increase  the  thickness 
of  the  membrane  to  0.7  mm.  The  external  membrane  was  bonded 
onto  a  specially  designed  membrane  clamping  ring.  This 
eliminated  one  source  of  leaks,  namely  through  the  0-rings 
clamped  onto  the  membrane  at  this  location.  The  thickness 
of  this  membrane  is  less  than  that  reported  in  Dusseault's 
paper.  This  membrane  is  approximately  2  mm  thick. 

A  detailed  description  of  sample  assembly  and  test 
procedures  is  given  in  a  later  section  of  this  chapter. 
Details  of  the  assembled  cell  and  measuring  systems  are 
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shown  on  Plate  2.1. 


2.4  Advantages  and  Disadvantages  of  the  HCTD 

The  hollow  cylinder  test  has  been  selected  for  this 
research  as  it  enables  control  of  all  three  principal 
stresses.  However,  this  is  also  possible  by  utilizing  a 
true  triaxial  cell.  Frydman( 1973 )  has  discussed  the 
problems  associated  with  this  apparatus:  end  effects  and 
strain  measurements.  The  end  effects  are  the  major  drawback 
of  the  true  triaxial  apparatus.  These  end  effects  are  not 
encountered  in  the  HCTD,  and  the  measurement  of  strain  is 
accomplished  by  measuring  the  volume  changes  occurring 
within  the  cell,  bore  and  sample. 

The  HCTD  has  additional  advantages  as  discussed  by 
Suklje  and  DrnovseM 1 965 ) : 

1.  The  stress  state  is  better  defined  without  uncontrolled 
stress  concentrations. 

2.  Tests  can  be  conducted  at  various  stress  states. 

3.  Deformation  anisotropy  can  be  taken  into  account. 

4.  Long-term  and  drained  tests  can  be  more  easily 
conducted . 

Despite  these  advantages,  there  still  exist  several 
drawbacks  to  the  apparatus: 

1.  A  non -homogeneous  stress  state  exists  in  the  sample 
which  complicates  the  interpretation  of  the  test 
results.  Assumptions  have  to  be  made  regarding  the 
stress  distribution  within  the  sample.  This  can  be  done 
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Plate  2.1  Instrumented  Cell  and  Measuring  System 
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by  using  elastic  theory  or  by  specifying  a  failure 
criterion  to  calculate  the  stresses  in  the  plastic  zone. 

2.  Special  rubber  membranes  are  required  to  resist  high 
cell  pressures  without  rupturing  or  leaking. 

3.  Problems  are  encountered  in  the  preparation  of  the 
samples  (discussed  in  more  detail  later). 

These  drawbacks  still  do  not  detract  from  the  versatility  of 
the  apparatus.  It  is  possible  to  conduct  plane-strain  tests 
along  different  stress  paths  in  a  relatively  simple  manner, 
and  with  accurate  control  of  pressures  and  volume  changes. 
Tests  may  be  run  either  in  an  undrained  mode  with 
pore-pressure  measurements,  or  with  full  drainage.  A  single 
test  yields  much  valuable  data: 

1.  Stress-strain  curves  based  on  average  stresses  within 
the  sample. 

2.  Volume  change  behaviour  due  to  dilatancy 

3.  Pore-pressure  response(for  undrained  tests). 

4.  Closure  of  the  bore  with  change  in  pressure. 

5.  Peak  strength  and  hence  strength  parameters. 

6.  Determination  of  both  elastic  parameters  (E  and  V  )• 

The  factor  which  contributed  largely  to  the  use  of  this 
apparatus  in  the  testing  program  is  the  fact  that  the  stress 
path  to  failure  in  a  shaft,  and  to  a  great  extent  in  a 
tunnel,  can  be  more  easily  and  closely  modelled  than  in  any 
other  piece  of  laboratory  equipment. 
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2.5  Outline  of  laboratory  program 

It  was  initially  intended  to  investigate  the  stress 
path  dependency  of  oil  sands  by  conducting  three  different 
series  of  drained  tests  in  the  HCTD.  The  laboratory  program 
was  modified  when  problems  with  the  HCTD  were  encountered 
(see  Appendix  F).  One  series  of  tests  were  conducted  in  the 
HCTD  (Series  A)  and  two  other  series  in  a  high  pressure 
triaxial  cell  (Series  B  and  C).  The  stress  path  followed  in 
each  of  the  three  series  of  tests  is  shown  in  Figure  2.1  and 
discussed  below. 

2.5.1  Hollow  Cylinder  Triaxial  Tests  (Series  A) 

This  series  consisted  of  two  plane-strain  drained  tests 
consolidated  to  7  Mpa  and  10  MPa  respectively.  The  samples 
were  sheared  by  maintaining  a  constant  external  pressure  and 
decreasing  the  internal  pressure  incrementally  until 
collapse  of  the  sample  occurred.  The  volume  changes  of  the 
internal  membrane  and  the  pore  fluid  were  monitored 
throughout  the  test. 

2.5.2  Triaxial  Compression  Tests  (Series  B) 

Two  samples  were  consolidated  to  7  MPa  and  10  MPa 
respectively  and  sheared  at  a  constant  strain  rate  with  full 
drainage  by  increasing  the  axial  stress  while  maintaining  a 
constant  cell  pressure.  Volume  changes  of  the  pore  fluid 
were  continuously  monitored  against  a  constant  back 


pressure . 
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Figure  2.1  Stress  paths  followed  in  the  laboratory  tests 
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2.5.3  Constant  Mean  Principle  Stress  Tests  (Series  C) 

In  these  tests,  two  samples  were  consolidated  to  5  MPa 
and  7  MPa  respectively.  The  samples  were  sheared  at  a 
constant  strain  rate  under  drained  conditions  with  constant 
mean  principal  stress  (  J1  =(  (J,  ♦  Gz  *  CJ3  )/3  =  constant  )  by 
increasing  the  axial  stress  and  decreasing  the  cell 
pressure.  Volume  changes  were  again  monitored  throughout 
the  test  against  a  constant  back  pressure. 

Details  of  the  sample  assembly  and  test  procedures  are 
given  in  a  later  section.  An  oedometer  test  was  used  to 
determine  the  compressibility  of  oil  sands  at  high  pressure. 

2.5.4  Selection  of  strain  rate 

The  strain  rate  selected  for  Series  B  and  C  tests  was 
based  initially  on  the  method  proposed  by  Bishop  and 
Henkel ( 1957 , pp 125- 127 )  where  the  time  to  failure  depends 
upon  both  the  drainage  conditions  and  the  coefficient  of 
consol idat ion . 

Negligible  volume  changes  are  occurring  in  the  sample 
after  10  minutes  for  each  application  of  cell  pressure.  The 
theoretical  t  jqq  value  is  less  than  10  minutes,  and  the 
calculated  time  to  failure  is  only  1  hour. 

The  coefficient  of  consol idat ion ,  Cv  ,  is  a  function  of 
both  the  coefficient  of  compressibility  and  the 
permeability.  The  coefficient  of  compressibility  is  stress 
dependent,  and  the  permeability  of  water  through  the  oil 
sand  sample  will  increase  as  the  sample  dilates. 
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The  value  selected  for  the  coefficient  of  consolidation 
was  obtained  from  the  isotropic  consolidation  test,  and  the 
calculated  time  to  failure  should  be  conservative.  However, 
due  to  the  variability  of  these  parameters  and  the  presence 
of  the  bitumen,  the  time  to  failure  in  these  tests  was  taken 
as  4  hours. 


2.6  Sampling  and  Sample  Preparation 

2.6.1  Sampling 

Sampling  of  oil  sands  is  difficult  as  sample 
disturbance  results  from  the  expansion  caused  by  gas  coming 
out  of  solution  when  the  confining  pressure  is  reduced. 
Dusseau 1 t ( 1 980 )  has  discussed  this  phenomena  in  detail  and 
has  shown  that  conventional  geotechnical  methods  of 
obtaining  samples  are  inadequate  for  use  in  the  oil  sands. 
The  best  methods  for  obtaining  relatively  high  quality 
samples  is  by  pressure  coring  and  in  situ  freezing 
procedures,  or  within  natural  slopes  where  the  gas  has  been 
depleted . 

In  situ  freezing  techniques  are  very  costly,  hence  it 
was  decided  to  obtain  block  samples  from  an  open  pit  mine 
where  the  gas  has  slowly  escaped  and  thus  disturbance  is 
minimized.  Preliminary  determination  of  sample  densities 
gave  a  bulk  density  of  2.0  g/cc  (porosity  of  35%)  which 
compares  favourably  with  the  results  reported  by 
Dusseau 1 t ( 1980 ) . 
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The  material  used  was  obtained  from  the  GCOS  (Great 
Canadian  Oil  Sands  Ltd.,  now  Suncor  Inc.)  open-pit  mining 
operation  north  of  Fort  McMur r ay .  The  samples  were 
collected  from  the  middle  McMurray  Formation  and  are 
bitumen-rich  oil  sands.  The  sampling  procedure  used  was  to 
cut  a  block  from  the  large  blocks  of  material  that  have 
fallen  off  the  face  of  the  excavated  walls  along  exfoliation 
fractures.  Immediately  after  a  block  had  been  cut,  it  was 
wrapped  in  a  plastic  bag  and  well  bound  with 
f ibre-rei nforced  tape.  The  sample  was  then  transported 
carefully  down  the  slope  and  stored  on  a  thick  foam  mattress 
in  the  back  of  the  vehicle.  Extreme  care  was  taken  in  the 
handling  of  the  samples  in  order  to  minimize  disturbance  and 
possible  breakage  of  the  larger  blocks. 

On  arrival  at  the  laboratory,  the  blocks  were  wrapped 
in  cheese-cloth  and  waxed  to  prevent  any  further  moisture 
loss  and  deterioration.  The  waxed  samples  were  then  stored 
in  a  moist  room  until  they  were  required  for  testing. 

2.6.2  Sample  Preparation 

2.6.2. 1  Oedometer  and  Tri axial  Samples 
All  triaxial  samples  were  obtained  from  a  single  block 
of  rich  oil  sand,  while  the  oedometer  sample  was  obtained 
from  a  block  of  lean  oil  sand.  The  block  was  removed  from 

O 

the  moist  room  and  allowed  to  freeze  in  a  cold  room  at  -18 
C.  The  block  was  cut  into  smaller  pieces  approximately  10 
cm  square  by  15  cm  long  with  a  diamond  saw.  These  blocks 
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were  wrapped  in  plastic  bags,  taped,  then  stored  in  the  cold 
room  in  a  container  of  dry  ice. 

When  a  sample  was  required  for  testing,  the  block  was 
unwrapped  and  placed  in  a  lathe  (installed  in  the  cold 
room).  The  sample  was  trimmed  using  tungsten-carbide  bits. 
It  is  important  that  the  sample  does  not  heat  up  during  this 
operation  as  this  will  result  in  some  expansion.  To 
minimize  this  effect,  the  sample  was  returned  several  times 
to  the  dry  ice  container  to  keep  it  well  frozen.  The  entire 
length  of  the  sample  is  never  trimmed  as  some  end  bearing 
surface  is  needed  to  clamp  the  sample  in  the  lathe. 

The  ends  of  the  trimmed  sample  are  cut  with  the  diamond 
saw  to  approximately  the  required  length  and  the  sample  is 
then  refrozen  in  the  dry  ice.  The  final  stage  of  the 
preparation  consisted  of  grinding  the  ends  of  the  sample  on 
a  belt  sander  to  the  correct  length.  A  metal  guide  was 
mounted  on  the  sander  to  ensure  that  the  ends  of  the  samples 
were  square  and  parallel.  Extreme  care  was  taken  to  ensure 
that  the  ends  did  not  heat  up. 

2. 6. 2. 2  Hollow  Cylinder  Samples 

The  preparation  of  these  samples  required  some 
experimentation.  The  method  described  below  is  that  which 
was  finally  used  in  the  preparation  of  the  actual  test 
specimens . 

A  block  sample  was  left  to  freeze  in  the  cold  room  and 
then  cut  into  a  cube  of  approximately  20  cm  square  by  25  cm. 
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This  cube  was  frozen  into  a  wooden  box  containing  saturated 
sand.  At  least  5  cm  of  sand  was  placed  beneath  the  sample 
and  about  5  cm  around  the  perimeter  of  the  sample.  The 
wooden  box  was  clamped  onto  a  milling  machine  table  (in  the 
cold  room)  and  dry  ice  packed  around  the  sample  to  its  full 
height . 

The  block  was  left  to  freeze  for  at  least  two  hours 
after  the  dry  ice  was  placed.  A  diamond  core  barrel,  51.4 
mm  external  diameter,  was  used  to  core  the  centre  hole. 
Coring  was  done  dry,  pressured  nitrogen  being  used  to  blow 
up  the  cuttings.  The  coring  was  done  slowly  and  no  more 
than  5  cm  cored  at  any  one  time  before  letting  any  built-up 
heat  dissipate.  Dry  ice  was  placed  in  a  bag  and  inserted  in 
the  cored  hole  to  assist  in  the  cooling  operation.  The  core 
barrel  had  to  be  constantly  cleaned  to  remove  the  broken 
core.  Inspection  of  the  bore  immediately  after  each  coring 
operation  revealed  no  signs  of  heating  on  this  inside  wall 
of  the  sample. 

A  similar  procedure  was  used  to  obtain  a  sample  of  15 
cm  external  diameter.  A  15  cm  diameter  core  barrel  was  used 
which  had  a  larger  clearance  than  most  core  barrels  on  the 
inside  wall.  This  facilitated  the  collection  and  removal  of 
cuttings.  The  barrel  was  centred  over  the  first  hole  and 
coring  carried  out  in  a  similar  manner.  The  trimming  of 
this  oversize  sample  was  done  on  the  lathe.  An  aluminum 
shaft  was  frozen  into  the  bore  and  filled  with  dry  ice 
before  mounting  it  in  the  lathe.  Trimming  was  carried  out 
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as  described  earlier  with  regular  refreezing  periods.  The 
ends  of  the  sample  could  be  finished  on  the  lathe  due  to  the 
mounting  procedure  used.  This  ensured  that  the  ends  were 
f lush  and  para  1  lei . 

The  sample  preparation  techniques  described  above  are 
illustrated  in  part  on  Plates  2.2  and  2.3. 


2.7  Sample  Assembly  and  Test  Procedures 

The  triaxial  samples  are  assembled  as  for  a  normal 
triaxial  test  on  any  other  material.  Once  the  sample  has 
been  placed  in  the  membrane,  the  membrane  is  filled  with 
water  to  displace  the  air.  The  sample  is  allowed  to  thaw 
against  a  cell  pressure  of  350  KPa  and  a  back  pressure  of 
200  kPa  prior  to  the  commencement  of  the  consolidation  test. 
The  consolidation  is  conducted  against  a  back  pressure  of 
1050  kPa,  each  increment  of  cell  pressure  being  maintained 
until  no  noticeable  volume  change  of  the  pore  fluid  is 
occurring  (usually  10  min.).  The  final  consolidation 
pressure  is  applied  for  at  least  30  min.  before  starting  the 
shear  test,  which  is  carried  out  at  a  constant  strain  rate. 
For  the  Jl-constant  tests,  a  programmable  calculator  was 
used  to  determine  the  cell  pressure  required  at  any 
particular  axial  stress.  The  cell  pressure  was  determined 
by  allowing  for  an  area  correction  for  the  axial  stress 
following  Bishop  and  Henke  1 ( 1 962 ) .  Throughout  the  test  the 
cell  pressure  was  calculated  at  1  to  2  minute  intervals  and 
the  pressure  reduced  as  required.  This  procedure  worked 
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Plate  2.2  Sample  preparation  techniques 
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Plate  2.3  Sample  preparation  techniques 
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well  as  the  maximum  variation  in  J1  was  less  than  2%. 

The  hollow  cylinder  samples  were  mounted  onto  the  base 
plate  in  the  cold  room.  The  external  membrane  clamp  is 
fastened  to  the  base  plate  and  the  sample  inserted.  The 
internal  membrane  is  inserted  and  the  ram  placed  on  top  of 
the  load  cap.  The  external  membrane  is  sealed  against  the 
load  cap  and  the  ram  by  means  of  0-rings  and  hose  clamps. 

The  base  plate  is  removed  from  the  cold  room  and  the  cell 
p 1 aced  onto  it. 

The  cell  is  filled  with  water  and  an  initial  pressure 
of  55  kPa  applied  to  both  membranes  prior  to  installing  the 
cell  in  the  plane-strain  frame  .  The  cell  pressure  is  raised 
slightly  while  the  system  is  flushed  with  a  pressure  of  70 
kPa  to  100  kPa.  This  operation  removes  air  from  the  porous 
stones  and  between  the  sample  and  membranes.  After 
flushing,  the  cell  and  bore  pressures  are  increased  to  350 
kPa  while  the  back  pressure  is  increased  to  200  kPa.  Under 
these  conditions  the  sample  is  left  to  thaw  overnight.  The 

consolidation  test  is  conducted  against  a  back  pressure  of 

/ 

1050  kPa.  The  shear  test  is  stress  controlled.  The  cell 
pressure  is  kept  constant  and  the  internal  pressure 
decreased  in  500  kPa  increments.  Each  increment  is 
maintained  for  at  least  5  min.  or  until  no  further  volume 
changes  occur.  The  increments  are  reduced  to  about  140  kPa 
in  the  final  stages  of  the  test. 
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2.8  Sample  Saturation 

The  method  of  testing  for  sample  saturation  as  outlined 
in  Bishop  and  Henkel (1962)  is  not  applicable  when  testing 
oil  sand.  The  low  compressibility  of  the  material  results 
in  a  B-value  less  than  unity  at  saturation.  The  samples 
tested  had  an  initial  degree  of  saturation  of  only  80%, 
hence  a  large  back  pressure  was  required  to  saturate  them. 

Full  saturation  can  be  accomplished  by  use  of  a  back 
pressure,  but  it  is  also  a  function  of  time  (Black  and 
Lee, 1973).  The  samples  had  an  applied  back  pressure  of  1050 
kPa  (150  psi )  and  this  pressure  should  be  sufficient  to 
fully  saturate  them.  After  application  of  the  back 
pressure,  the  pore  volume  change  was  monitored  and  no 
noticeable  volume  changes  were  occurring  prior  to  the  start 
of  the  consolidation  test.  However,  the  sensitivity  of  the 
burettes  may  not  have  been  sufficient  for  this  purpose.  The 
time  interval  from  the  application  of  the  back  pressure  to 
the  start  of  the  shear  test  was  usually  1.5  hours,  and  this 
should  also  assist  in  saturating  the  samples. 

Full  saturation  is  critical  when  conducting  undrained 
tests  and  a  satisfactory  laboratory  procedure  to  test  for 
saturation  will  be  required.  Undrained  tests  on  oil  sand 
will  be  attempted  by  another  researcher  soon,  and  two 
methods  of  testing  for  saturation  are  outlined  below. 

The  first  method  is  based  on  the  fact  that  the  B-value 
is  independent  of  back  pressure  ( Wi ssa , 1 969 ) .  The  back 
pressure  is  increased  for  each  B-test  while  keeping  the 
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effective  consolidation  pressure  constant.  Full  saturation 
results  when  the  calculated  B-value  remains  constant.  In 
order  to  maintain  a  constant  effective  consolidation 
pressure,  each  increment  of  back  pressure  and  cell  pressure 
must  be  exactly  equal  before  measuring  B.  Another  method 
proposed  by  Wissa(1969)  is  to  apply  several  large  increments 
of  cell  pressure  without  allowing  drainage.  The  B-values 
are  measured  and  will  be  the  same  (or  slightly  smaller)  for 
each  increment  of  cell  pressure  at  saturation. 

Other  methods  may  be  devised,  but  the  method  chosen 
should  be  relatively  simple  and  quick.  This  procedure  may 
well  be  the  critical  part  of  any  undrained  tests  that  are 
performed,  and  serious  attention  must  be  directed  towards 
this  problem. 


2.9  Measurements  and  Measuring  Systems 

All  data  collected  during  the  test  were  recorded  by 
electronic  measuring  systems  except  for  the  volume  changes. 
Axial  loads  are  measured  on  a  89  kN  load  cell  and  axial 
strain  with  a  linearly  variable  displacement  transducer 
(LVDT).  The  cell  pressures  and  back  pressures  are  measured 
with  13.8  MPa  pressure  transducers.  Measurements  of  volume 
changes  are  made  on  reversible  burette  type  volume  change 
indicators.  For  the  hollow  cylinder  tests  these  indicators 
have  an  accuracy  of  0.1  cc.  The  triaxial  samples  require  a 
more  sensitive  burette  due  to  the  smaller  volume  of  these 
samples,  and  burettes  with  an  accuracy  of  0.01  cc  were  used. 
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The  pressure  systems  used  in  the  laboratory  are 
discussed  in  Appendix  F,  which  also  deals  with  the  problems 
encountered  in  the  laboratory  during  this  research  program. 


CHAPTER  3 


CALCULATION  OF  STRESSES  AND  STRAINS  IN  THE  HOLLOW  CYLINDER 
3.1  Introduction 

The  stress  distribution  within  the  hollow  cylinder 
depends  upon  both  the  stress  state  and  the  type  of  test 
being  analysed.  If  the  state  of  stress  is  such  that  no 
yielding  occurs  at  any  point  within  the  cylinder,  then  the 
stress  distribution  can  be  determined  from  the  theory  of 
elasticity.  When  yielding  occurs,  plasticity  solutions  are 
required  to  determine  the  distributions. 

The  series  of  tests  conducted  on  the  hollow  cylinder 
samples  were  performed  under  plane-strain  conditions.  For 
this  case,  the  axial  stress  is  the  intermediate  principal 
stress  and  the  sample  deforms  axi symmetr ica 1 ly ,  the  slip 
lines  occurring  in  the  r-0  plane.  The  distribution  of 
stresses  can  be  determined  from  limiting  equilibrium 
conditions  by  specifying  a  failure  criterion  for  the 
material.  The  solution  for  the  case  where  the  axial  stress 
is  the  major  principal  stress  requires  that  a  stress-strain 
law  be  specified  to  solve  for  the  stress  distribution  ( Wu  et 
al . , 1963) . 

In  this  series  of  tests,  the  external  pressure  is  Kept 
constant  and  the  internal  pressure  decreased  until  the 
sample  collapses.  The  tangential  stress  is  the  major 
principal  stress  and  the  radial  stress  the  minor  one.  The 
stress  state  is  such  that  yielding  first  occurs  at  the 
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internal  wall  of  the  sample. 

The  solutions  for  the  distribution  of  radial  and 
tangential  stresses  within  the  hollow  cylinder  are  given  in 
the  following  sections  for  both  the  case  of  no  yielding  and 
for  the  case  when  yielding  has  occured. 


3.2  Elastic  Solution  for  the  Stress  Distribution 

The  axi symmetric  problem  may  be  solved  in  terms  of  the 
equilibrium  equations,  Hooke's  Law,  the  compatibility 
requirements  and  the  Airy  stress  function.  By  applying  the 
appropiate  boundary  conditions,  the  solution  for  the  radial 
and  tangential  stress  distributions  are  determined  : 


...3.1(a) 


(Je=-^r  +  2C 


_ 3.1(b) 


where : 


....3.1(c) 


The  symbols  used  in  the  above  equations  are  illustrated  in 
F igure  3.1. 
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a  <  r  <  b 


Figure  3.1  Definition  of  the  terms  used  in  the  Stress 
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The  derivation  of  these  equations  is  given  in  Appendix 
A.  The  average  stresses  in  the  sample  are  found  by 
integrating  the  above  stress  distributions  across  the  width 
of  the  cylinder  and  dividing  by  the  wall  thickness.  The 
results  are  given  below  : 

( CT  )  s  ?r  +  A  ....3.! 


( ^0  )ave  58  2C  - 
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ab 


....3.3 


3.3  Stress  Distribution  after  Yielding 

As  discussed  before,  yielding  commences  at  the  internal 
wall  of  the  sample.  A  plastic  zone  develops 
axi symmetrical ly  around  the  inner  surface,  while  the 
remainder  of  the  sample  is  still  behaving  elastically.  This 
plastic  zone  is  shown  di agramat ical ly  in  Figure  3.1.  The 
stresses  in  the  elastic  zone  cannot  be  calculated  from  the 
equations  given  previously  as  the  boundary  conditions  have 
changed.  The  solution  for  the  stresses  in  the  elastic  zone 
requires  that  the  extent  of  the  plastic  zone  be  known. 

3.3.1  Stresses  in  the  Plastic  Zone 

The  Mohr -Coulomb  failure  criterion  is  assumed  to  be 
valid  for  this  material.  The  validity  of  this  assumption  is 
discussed  in  Chapter  4.  By  combining  the  equilibrium 
equations  with  the  failure  criterion,  the  stress 
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distribution  is  obtained  : 


CTr  =  J^[ArN'’  -  2cN1/2  ] 
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where  A  is  solved  from  the  boundary  conditions  : 


A  =  — 


P;  (N-1 )  +  2  c  N 


1/2 


N-1 


N  - 


1+sin  <p 

1  -  sin  0 


The  derivation  of  these  equations  is  given  in  Appendix  B. 

3.3.2  Stresses  in  the  Elastic  Zone 

The  solution  for  these  stresses  can  be  obtained  from 
Equations  3.1  provided  that  the  boundaries  of  the  elastic 
zone  are  modified.  The  elastic  zone  is  now  defined  by  an 
internal  radius  (R)  and  an  external  radius  (b).  The 
pressures  on  the  internal  and  external  walls  are  Op  and  p0 
respectively. 

Substitution  of  these  values  into  Equation  3.1  gives  : 


cr; 


k  ~  Op  R 
fcf-R* 


RVtP.-CTp) 

rNtf-R1) 


....3.6(a) 
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-  Op  R* 

rfr-  * 


.  . . .3.6(b) 
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It  is  necessary  to  calculate  the  radius  of  the  plastic 
zone  and  the  radial  stress  at  the  elastic~plastic  boundary 
before  the  stress  distributions  can  be  determined. 

•  3.3.2. 1  Radius  of  the  Plastic  Zone 

The  radius  is  determined  by  assuming  that  the  radial 
stress  at  the  boundary  is  continuous.  Before  yielding 
occurs,  the  sum  of  the  elastic  radial  and  tangential 
stresses  (  CT0  )  is  a  constant  value  at  any  point  within  the 
cylinder  (for  each  set  of  applied  boundary  pressures). 

After  yielding  has  occurred,  the  sum  of  the  elastic  stresses 
in  the  elastic  zone  is  also  a  constant  value,  but  this 
constant  is  a  function  of  the  radius  of  the  plastic  zone. 

The  solution  for  the  radius  of  the  plastic  zone  is 
simplified  if  it  is  assumed  that  this  constant  is  of  the 
same  magnitude  as  would  be  calculated  if  no  yielding  had 
occurred.  In  Chapter  6,  the  validity  of  this  assumption  is 
discussed . 

The  radius  of  the  plastic  zone,  based  on  the 
assumptions  outlined  above,  is  determined  by  equating  the 
values  of  the  radial  stress  at  the  boundary  calculated  from 
the  elastic  and  plastic  solutions.  The  solution  is  given 
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below  and  the  derivation  shown  in  Appendix  B. 


R 


a 


CT0(N-n  +  2cN1/2 

Pj(N-1  )  +  2cN1/2 


«  •  •  •  W  > 


The  radial  stress  at  the  elastic-plastic  boundary  is 
found  from  Equation  3.4  : 


3.4  Average  Stresses  within  the  Hollow  Cylinder 
3.4.1  Average  Tangential  Stress 
The  average  tangential  stress  can  be  found  by 
integrating  the  stress  distribution  across  the  width  of  the 
wall  and  dividing  by  the  wall  thickness.  However,  the 
average  stress  can  also  be  found  by  considering  the  static 
equilibrium  of  a  half-section  of  the  cylinder.  Figure  3.2 
shows  this  section  and  the  calculation.  This  approach  is 
simpler  and  the  solution  is  independent  of  the  extent  of  the 
plastic  zone.  The  average  tangential  stress  is  : 


l  t  I  •  V  ( 


3.4.2  Average  Radial  Stress 

The  average  radial  stress  can  be  determined  in  two 
ways.  A  linear  distribution  of  the  stresses  may  be  assumed 
(or  any  other  shape)  and  the  average  stress  calculated  from 
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Ft 


For  equilibrium  : 


2Ft  +  p.  (2a)  -  Pe(  2  b)  -  0 


R  = 


=  Peb  - 


PiQ 


<cu 


e'ave 


b-a 


Pgb  -  P;Q 

b-a 


Figure  3.2  Static  Equilibrium  conditions  for  the  half 

Cylinder 
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distribution,  or  the  stress  distribution  given  in 
Equations  3.4  and  3.6(a)  may  be  integrated  across  the  width 
of  the  cylinder  and  divided  by  the  wall  thickness. 

The  average  radial  stress  based  on  a  linear 
distribution  is  : 


( lave 


Pfi±Pi 

2 


3.10 


The  average  stress  found  by  integration  is  derived  in 
Appendix  B  and  given  for  a  cohesion  less  material  by  : 


\ive 


_Ejl 

b  -  a 


....3.11 


where : 


Fr 


RN 

NaN-1 


a_  (b-R)  RN  1 
N  (  b+  R  )  a  N-1  J 


Peb 


(b-R) 

(b+R) 


Both  of  these  solutions  have  been  used  in  the  analysis 
of  the  hollow  cylinder  test  data.  Their  effect  on  the 
stress-strain  curves  and  strength  parameters  are  discussed 
in  Chapter  6. 


3.4.3  Average  Axial  Stress 

The  load  applied  at  the  top  of  the  sample  is  measured 
throughout  the  test  and  thus  the  average  axial  stress  can  be 
determined  directly  : 


V 


easirS  f«ixA  aesievA  £ .A. 
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(CTZW 


P  -  P;  A; 


_ 3.12 


The  meaning  of  the  symbols  are  as  defined  in  the 
nomenclature . 


3.5  Determination  of  the  Strength  Parameters 

The  strength  of  a  cohesionless  material  at  high 
confining  pressures  may  be  defined  by  a  secant  value  of  the 
friction  angle,  •  This  angle  is  determined  from  the 
effective  principal  stresses  at  failure  (for  a  Mohr-Coulomb 
mater i al ) : 


Oi'-  Oa 

CT/+  (J31 


3.13 


The  value  of  0S  determined  by  this  method  will  be 
dependent  on  the  magnitude  of  the  average  radial  stress,  and 
the  calculation  of  this  stress  will  have  a  large  influence 
on  the  friction  angle.  An  alternate  method  of  determining 
the  peak  friction  angle  is  found  from  the  elasto-plastic 
analysis.  At  collapse  of  the  sample,  the  radius  of  the 
plastic  zone  becomes  equal  to  the  external  radius  (b)  and 
the  radial  stress  will  thus  be  equal  to  the  external 
pressure  (  Pe  )  in  Equation  3.8.  Rearrangement  of  this 
equation  under  these  conditions  give  : 


I ,  j 
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.3.14 


The  pressure  ratio  and  ratio  of  the  radii  at  collapse 
are  Known  from  the  test  results,  and  the  strength  parameter 
(N)  can  be  calculated. 

Comparison  of  the  friction  angles  determined  by  these 
methods  are  given  in  Chapter  6. 


3.6  Calculation  of  the  Average  Strains  in  the  Cylinder 

The  calculation  of  the  strains  are  based  on  average 
strain  increments  rather  than  total  strains.  Strain 
increments  have  been  used  as  the  test  is  performed  by 
reducing  the  internal  pressure  incrementally.  Total  strains 
are  calculated  by  adding  the  strain  increments. 

The  principal  strain  axes  are  assumed  to  coincide  with 
the  principal  stress  axes,  and  compressive  strains  are 
considered  as  being  positive.  The  small  strains  theory  is 
used ,  i . e  : 


dv  -  d€z  +  d£r  +  d€*  - J  -  1 

=  d€r  +  d£*  for  plane-strain 

The  calculation  of  strains  by  this  method  has  been 
dealt  with  by  Procter ( 1 967 ) .  His  derivations  are  presented 
in  a  form  that  can  be  readily  programmed.  A  copy  of  these 
derivations  is  given  in  Appendix  D. 
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3.6.1  Average  Radial  Strain 

The  average  strain  increment  is  defined  as  : 


d€ 


r 


mean  change  in  specimen  wall  thickness 

mean  wall  thickness 


The  wall  thickness  at  any  increment  is  calculated  from  the 
measured  volume  changes  of  the  bore  and  sample  : 


The  average  radial  strain  increment  is  : 


d€ 


r 


(b;-bj)  -  (a;  gj)1  x  2 
(b;  +  bj)-  (Qj+aj)J 


.3.16 


The  jth  increment  is  the  current  increment  and  the  ith 
increment  the  previous  one.  The  total  radial  strain  is 
found  by  summation  of  the  strain  increments  : 

(€r)j  *  £  d£r  ...  .3. 

J  1 

3.6.2  Average  Volumetric  Strain 

The  average  volumetric  strain  increment  is  given  by  : 


'  ■ 
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(AvS)j  -  (A  Vs); _ 1  ^  2  _ 3.18 

2  ( Vs  )0  -  ( AVs)j  -(DVsljJ  *  2 


The  total  volumetric  strain  is  : 

(v)j  =^dv  ....3.19 

1 

3.6.3  Average  Tangential  Strain 

When  the  average  radial  and  volumetric  strain 
increments  have  been  calculated,  the  average  tangential 
strain  increment  is  determined  from  Equation  3.15  : 

d£0  8  dv  -  d€r  ....3.20 


The  total  tangential  strain  is  : 

(6*)j-id€*  •••■3.21 

1 

3.7  Determination  of  the  Elastic  Deformation  Moduli 

Young's  Modulus  (E)  is  obtained  from  the  linear  elastic 
constitutive  equations  : 


A €r  -  Act.  +  Acrz)] 


•  *  • 


3.22(a) 
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Ae,  =y[Ao; -v{&(jr  +Acr2)]  ••••3.22(b) 

A«z  -  Aoi  -WAcTr  +  Aa*)]  •"•3'22(c) 


For  plane-strain  conditions,  Equation  3.22(c)  is  set  equal 
to  zero  and  the  following  relationship  is  obtained  : 


A  o'*  =  +A  cr9) 


Poisson's  ratio  is  determined  from  Equation  3.23  : 

,,  _  A  Gz  _ 3.24 

Acrr+Acr. 

Young's  Modulus  is  calculated  by  substituting  Equation  3.23 
into  Equation  3.22(b)  : 

E  ACT.-(T^)AO~r  . . . .3.25(a) 

Eps  “  1  -i>r  ' 

.  A(Jd 

A^e 


a  a  a  a 


3.25(b) 
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In  the  hollow  cylinder  test,  yielding  occurs  at  the 

p 

inside  wall  at  a  pressure  ratio  (-=£-)  of  only  2  to  3.  The 

”i 

elastic  modulus  is  thus  defined  for  only  a  small  portion  of 
the  average  stress-strain  curve.  The  use  of  a 
psuedo-elastic  modulus  (defined  as  the  slope  of  the  average 
stress-strain  curve  at  any  stress  level)  is  dealt  with  in 
Chapter  4. 


3.8  Summary 

The  calculation  of  stresses  and  strains  within  the 
hollow  cylinder  has  been  presented.  An  incremental  approach 
has  been  used  for  the  calculations  of  strains  due  to  the 
manner  in  which  this  series  of  tests  have  been  conducted. 
Several  assumptions  have  been  made  in  the  analysis  of  the 
stress  distributions.  Some  of  these  assumptions  are  those 
typically  made  for  soils,  the  others  for  simplicity  of  the 
solution  alone.  In  either  case,  the  validity  or  influence 
of  the  assumptions  on  the  analysis  of  the  test  data  will  be 
discussed  in  the  following  chapters. 


« 
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CHAPTER  4 


DEVELOPMENT  OF  THE  NUMERICAL  MODEL 
4.1  Introduction 

The  interpretation  of  data  from  a  hollow  cylinder 
triaxial  test  is  a  problem  that  has  as  yet  not  been  properly 
examined.  Assumptions  have  to  be  made  regarding  the  stress 
distribution  within  the  cylinder.  These  assumptions  have 
been  discussed  in  Chapter  3. 

Once  a  stress  distribution  has  been  determined,  the 
problem  of  how  to  interpret  the  elastic  deformation  moduli 
remains.  Linear  elastic  formulations  enable  the  Young's 
modulus  to  be  defined  in  terms  of  deviatoric  stress  and 
tangential  strain  as  presented  in  Chapter  3.  It  still 
remains  to  be  shown  that  the  moduli  determined  in  such  a 
manner  from  the  hollow  cylinder  data  do  in  fact  represent 
the  true  moduli  for  the  material. 

The  writer  has  approached  this  problem  by  determining 
the  elastic  constants  (assuming  a  linear  radial  stress 
distribution)  from  the  average  stress-strai n  curve  for  the 
test.  These  moduli  are  then  used  as  input  in  an  idealized 
model  in  an  attempt  to  predict  the  deformation  of  the 
cylinder  that  are  measured  during  a  test.  The  analytical 
model  developed  for  this  problem  is  dealt  with  in  the 
following  sections.  The  stress  and  strain  determinations 
from  the  hollow  cylinder  test  have  been  described  in  the 
previous  chapter. 


46 


•  ", 


. 


47 


4.2  Formulation  of  the  Problem 

Oil  sands,  as  well  as  other  materials,  exhibit 
curvilinear  failure  envelopes.  This  initiated  the 
development  of  an  axi -symmetr i c  solution  that  could 
incorporate  this  type  of  failure  envelope.  Curvilinear 
failure  envelopes  have  not  been  used  in  any  of  the 
elasto-plastic  closed  form  solutions  available  in  the 
1 i terature . 

The  strain-weakening  behaviour  of  oil  sands 
( Dusseau 1 t , 1 977 )  was  another  design  factor  to  be  considered. 
The  solution  should  be  able  to  model  such  behaviour  as 
closely  as  possible.  It  has  been  shown  that 
strain-weakening  behaviour  can  be  modelled  with  an 
e 1 ast i c-per feet ly  plastic  model  ( Korbi n , 1 976 ) . 

The  effect  of  confinement  on  the  Young's  modulus  is 
evident  from  the  test  results  which  show  a  marked  non-linear 
stress-strain  curve  (see  Chapter  5).  This  non-linear 
stress-strain  curve  is  the  result  of  the  influence  of 
reduced  confinement  of  the  sample  as  well  as  the  non-linear 
behaviour  of  the  material.  It  is  thus  necessary  to 
incorporate  the  strain-softening  effect  of  the  deformation 
modulus  into  the  model. 


4.3  Theoretical  Aspects  of  the  Model 

The  use  of  a  linear  or  non-linear  elastic  model  will 
not  be  satisfactory  in  the  prediction  of  deformations. 
These  models  work  well  at  small  strains,  but  in  a  hollow 


' 
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cylinder,  the  strains  are  large,  especially  once  yielding 
has  occurred. 

Elastic  models  have  further  limitations  (Hoeg,1978). 
They  do  not  model  typical  soil  behaviour  such  as  dilation 
(for  dense  sands),  stress  history,  stress  path  dependency 
and  s t r a i n -weaken i ng .  These  effects  can  be  modelled  by 
using  elasto-plastic  models. 

The  theory  of  plasticity  was  originally  developed  for 
metals  and  has  been  adapted  and  extended  to  model  soil 
behaviour.  The  deformations  are  determined  by  defining  a 
flow  rule  and  a  work-hardening  law,  which  are  relationships 
between  the  plastic  strain  increment  and  the  stress 
increment  vectors. 

The  undrained  behaviour  of  soils  can  still  be  modelled 
on  the  basis  of  classical  plasticity,  but  the  volume  changes 
that  occur  during  drained  behaviour  need  to  be  incorporated 
into  the  theory.  The  Mohr-Coulomb  strength  relationship  can 
be  considered  to  be  a  yield  criterion  and  a  flow  rule  must 
be  postulated  which  will  predict  the  plastic  volume  changes 
accurately. 

The  use  of  an  associated  flow  rule  with  the 
Mohr-Coulomb  yield  criterion  has  been  found  to  give 
unrealistically  large  plastic  volume  changes  (Hoeg,1978  ; 
Morgenstern , 1 975 )  and  this  led  to  the  introduction  of  a 
non-associ ated  flow  rule:  i.e.,  the  yield  criterion  and  the 
plastic  potential  function  are  not  identical.  The 
non-associ ated  flow  rule  can  account  for  volume  changes  that 
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occur  between  two  limits,  that  of  no  volume  change  and  that 
of  maximum  volume  change  (associated  flow  rule).  The  model 
developed  below  has  the  capability  of  incorporating  a 
dilation  parameter  (which  is  constant  within  the  plastic 
zone)  that  has  some  value  between  these  two  limits. 


4.3.1  Failure  Criterion 

The  Mohr -Cou lomb  failure  criterion  has  been  shown  to  be 
applicable  for  drained  tests  on  dense  sand  (Kirkpatrick, 
1957).  -The  non-linearity  of  the  Mohr  envelope  is  modelled 
by  using  an  empirical  criterion  postulated  by 
Dusseaul t ( 1 977 )  for  cohesionless  materials  : 


T 


_ 4.1 


where  CTn  is  the  normal  stress  on  the  failure  plane 
and  a,  b  are  curve  fitting  parameters 
For  simplicity,  it  is  assumed  that  this  equation 
describes  the  failure  criterion  which  passes  through  the  top 
of  the  Mohr  circles,  i.e.  that  failure  occurs  along  planes 
of  maximum  shear  stress  (Figure  4.1).  The  normal  stress  is 
then  given  by  : 


On  = 


CTe+Or 


_ 4.2 


■ 


SHEAR  STRESS 


50 


Figure  4.1  Failure  criterion  adopted  in  the  model 


Figure  4.2  Stress  path  followed  by  a  typical  soil  element  in 

the  model 
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In  the  hollow  cylinder  test,  plane  strain  is  maintained 
in  the  vertical  direction  (in  this  series  of  tests).  The 
author  believes  that,  as  a  result  of  this  constraint, 
failure  planes  are  unlikely  to  be  the  same  as  for  the 
conventional  triaxial  test,  and  that  failure  along  planes  of 
maximum  shear  stress  is  most  probable.  However,  as 
mentioned,  it  is  simple  to  incorporate  any  other  analytic 
criterion  into  the  model.  The  Mohr-Coulomb  equation  could 
be  modified  to  give  the  following  criterion  : 

<7,  =  a  CTrb 


The  maximum  shear  stress  is  given  by  : 


Tn 


max 


2 


and  the  failure  criterion  may  be  rewritten  as  : 


^Tr  ( &r  \  ....4.3 

2  '  2  ' 

4.3.2  Determination  of  Stresses  in  the  Cylinder 

In  order  that  the  failure  criterion  is  not  violated, 
the  stress  path  followed  by  an  element  in  the  cylinder  is  as 
illustrated  in  Figure  4.2.  The  stress  path  AB  represents 
the  stress  increments  while  the  cylinder  is  behaving 
elastically,  and  Point  B  represents  the  stress  state  at 
yielding  of  the  element. 

When  the  element  has  yielded,  the  stresses  are  defined 
by  the  failure  criterion  and  hence  decrease  with  decreasing 
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internal  pressure  (path  BC).  The  calculation  of  the 
stresses  is  dealt  with  in  the  following  section. 

4.3.2. 1  Pre- Yield  Condition 

The  stresses  in  the  cylinder  are  calculated  from  the 
linear  elastic  equations  derived  for  a  hollow  cylinder  (see 
Chapter  3).  The  stress  distributions  are  given  by  Equations 
3.1(a)  to  3.1(c). 

4. 3. 2. 2  Post -Yield  Condition 

The  stresses  in  the  plastic  zone  are  calculated  by 
considering  the  equilibrium  equation  and  failure  criterion 
in  the  plastic  zone  (Equations  B.2  and  4.3).  The  solution 
of  these  equations  is  complex  and  makes  a  closed  form 
solution  unsuitable.  A  finite-difference  technique  is  used 
to  solve  the  equations. 

For  the  axi -symmetric  case,  (Tr  =  fct(r),  and  the 
Taylor  series  for  this  function  can  be  expanded  as  follows  : 


2 


d  Or 


dr 


Ignoring  second  order  and  higher  terms  yield  : 


d  Or  _  Oq  Or 


From  Eqn.  B.2  : 


dr 
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The  radius  of  the  plastic  zone  is  calculated  by 
assuming  continuity  of  the  radial  stresses  at  the 
elastic-plastic  boundary.  At  the  boundary ,  CJr  =  (Jp  and  the 
failure  criterion  can  also  be  applied. 

From  Equations  3.6(a)  and  3.6(b)  we  get  : 


. . . .4.6(a) 


_ 4.6(b) 


Substituting  Equation  4.6(b)  into  the  failure  criterion 


gives  : 


2 


Gx  *  Oo  -  0  G0 


and : 


4  7 
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Equations  4.5  and  4.7  are  solved  iteratively.  The 
cylinder  is  divided  into  a  number  of  discrete  nodes.  The 
elastic  stresses  on  the  inside  wall  are  calculated  for  each 
load  increment  until  the  failure  criterion  must  be  invoked. 
At  this  stage,  the  tangential  stresses  are  related  to  the 
radial  stresses  by  the  failure  criterion.  These  stresses  on 
the  wall  form  the  first  set  of  values  to  be  used  in  the 
iterative  process .  ^  (Jr  is  calculated  for  the  increment  , 
and  thus  the  value  of  ((Jr).  ^  is  determined.  ( 0^).  ^  is 
found  from  Equation  4.3. 

The  value  of  ((Jr).  ^  and  (r)j^  are  substituted  into 
Equation  4.7  and  the  value  of  (7re  is  compared  to  CJp  .  The 
procedure  continues  until  either  the  difference 

j((Jp  -  CTre)j  is  within  the  specified  tolerance  or  the 
difference  ((Jp“  CTre)  changes  sign  from  one  iteration  to  the 
next.  At  this  point,  (Tr  and  r  correspond  to  the  radial 
stress  at  the  plastic  boundary  and  the  radius  of  the  plastic 
boundary  respectively. 

The  stresses  in  the  elastic  zone  are  then  calculated 
from  Equations  3.6(a)  and  3.6(b). 

4.3.3  Strain  Distribution 

4.3.3. 1  Elastic  Strains 

Before  yielding  occurs,  the  strains  are  determined  from 
Hooke's  law  and  the  elastic  stress  distribution  (see 

Equations  3.22(a)  to  3.22(c)). 

Strains  are  dealt  with  in  an  incremental  manner.  At 


r  :ss»  -10*  ? u.  I'Uw  ebranr  *4J  r^o  aaraa^Js  -?*■??* 
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the  completion  of  each  load  increment,  the  model  boundaries 
are  modified  by  the  amount  of  radial  displacement  at  each 
wall.  Integration  of  the  strain-displacement  equations 
yield  : 


**  ur  =  r^o 

ie'  Aur=rA€,  ....4.8 

This  incremental  strain  procedure  was  instituted  in  the 
model  since  the  loading  is  applied  incrementally,  and  the 
stress  distribution  is  sensitive  to  the  values  of  the 
boundary  radii . 


4. 3. 3. 2  Strains  in  the  Plastic  Zone 

The  strains  in  the  plastic  zone  are  assumed  to  consist 
of  an  elastic  and  plastic  portion,  thus  : 

s  .  .  .  .4.9(a) 


% . . .4.9(b) 


The  elastic  portion  of  the  strain  is  assumed  to  obey 
Hooke's  law  while  the  plastic  strains  are  defined  by  a  flow 
rule.  Both  the  radial  and  tangential  elastic  strains 
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(6 ®  i  ^r)  are  taken  as  being  constant  in  the  plastic  zone 
and  equal  in  magnitude  to  the  strains  at  the  elastic-plastic 
boundary,  which  simplifies  the  solution.  The  elastic 
strains  in  the  plastic  zone  can  thus  be  defined  by  : 


i-V 

E 


[Acrp  (1_^)AO«r] 


_ 4.10(b) 


The  strain  compatibility  equation  may  be  written  : 


r£e  +€P-6P  -€t- 


d 


3. 


dr 


_ 4.11(a) 


d£* 

Since  6®  is  constant  in  the  plastic  zone,  9  »0  > 


dr 


and  : 


d€| 


r  — — ~  —  £r  * 

dr  9 


_ 4.11(b) 


The  flow  rule  adopted  is  a  modified  form  of  the 
associated  flow  rule  ( Guenot , 1 978 )  and  is  given  as  . 


€pr  +<*€?  =  0 


_ 4.12 
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where  o C  =  dilation  parameter.  When  cK.  =1,  the  condition  of 
no  volume  change  is  modelled.  The  associated  flow  rule  is 
obtained  when  c<  -  N,  the  shear  strength  parameter.  This 
represents  the  case  of  maximum  volume  change. 

The  compatibility  equation  is  rewritten  as  : 

deg  \e*  -el  -€,p(i+°c)l  4  i3 

dr  r 


This  equation  can  also  be  solved  by  a  finite-difference 
technique,  similar  to  that* given  previously  for  the  stress 
distribution  : 

A<s»  .  (eSL  - (€S)i 

Ar  "  Ar 


A€e  ’  [(€r  -  6j)j  -  (€j).(l+«x)] 


•  •  •  • 


4.  14 


The  solution  is  obtained  by  selecting  an  arbitrary 
value  for  (c£)  at  the  inside  wall.  The  iteration  is 
carried  out  until  r  =  R  .  At  this  point,  the  plastic  strain 
must  be  zero.  The  value  of  (€«)•  is  changed  until  the 

solution  satisfies  the  tolerance  specified. 

% 

The  total  strains  are  the  sum  of  the  elastic  and 
plastic  strains.  The  plastic  strains  calculated  to  satisfy 
the  compatibility  equation  represent  the  total  plastic 
strain  at  any  stage  of  the  test.  These  values  are  then 
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transformed  into  plastic  strain  increments  by  subtracting 
the  total  plastic  strain  at  the  previous  load  increment. 

The  change  in  radial  displacement  is  calculated  from  the 
strain-displacement  equation  and  the  model  boundaries  are 
adjusted  as  described  earlier. 

4.3.4  Strain-Softening  Effect 

The  incremental  procedure  used  in  the  analysis  is  well 
suited  for  incorporating  a  deformation  modulus  that  is 
stress-dependent.  It  remains  to  be  ascertained  which 
particular  function  best  describes  this  softening  effect. 

An  average  stress-strai n  curve  from  a  hollow  cylinder 
test  on  Ottawa  sand  is  shown  in  Figure  4.3,  based  on  the 
assumption  of  a  linear  radial  stress  distribution.  It  can 
be  seen  from  this  curve  that  the  behaviour  of  the  material 
is  non-linear,  and  thus  the  tangent  Young's  modulus  (E)  will 
decrease  with  increasing  deviatoric  stress.  In  other  words, 
the  modulus  will  decrease  with  a  decrease  in  the  average 
radial  stress. 

A  suitable  function  has  to  be  defined  which  can  best 
describe  the  variation  of  the  deformation  modulus  within  the 
elastic  region  of  the  cylinder.  Three  alternate  functions 
were  examined  : 

1.  A  variation  of  the  modulus  with  average  radial  stress 

2.  A  variation  of  the  modulus  with  average  deviatoric 
stress 

3.  A  variation  of  the  modulus  with  the  radial  stress  at  the 
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Figure  4.3  Average  stress-strain  curve  for  the  Ottawa  sand 

test 
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elastic-plastic  boundary 

The  plots  of  these  functions  for  the  Ottawa  sand  test 
are  shown  in  Figure  4.4.  The  values  of  E  are  calculated 
from  the  slope  of  the  stress-strain  curve  (Figure  4.3)  at 
different  stress  levels.  The  wall  closures  are  compared 
with  the  actual  closures  in  Figure  4.5. 

The  function  that  yields  the  best  result  is  the  first 
one  (E  varies  with  average  radial  stress).  This  function 
can  be  described  by  a  bi-linear  distribution  of  the 
following  form  : 


....4.15(a) 


where  (Jpc  =  average  radial  stress  at  the  intersection  of 
the  above  lines. 

The  average  radial  stress  used  in  Equations  4.15 
depends  on  whether  yielding  has  occurred  or  not.  Prior  to 
yielding,  the  average  radial  stress  is  defined  as  : 


P; 


....4.16 


. 
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Fiqure  4.4  Possible  functions  for  defining  E  in  the  elastic 

zone 
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Figure  4.5  Influence  of  the  E-function  on  the  wall  closure 
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Once  yielding  has  occurred,  the  average  radial  stress 
is  given  by  : 


_ 4.17 


This  value  is  used  as  the  elastic  strains  in  the 
plastic  zone  are  determined  from  the  stresses  at  the 
elastic-plastic  boundary.  The  modulus  used  reflects  the 
average  confinement  in  the  elastic  zone. 

4.4  Output  from  the  Program 

The  model  calculates  and  outputs  the  following  data  : 

-  initial  input  parameters 

-  radial  and  tangential  stress  at  each  node 

-  value  of  the  radius  at  each  node 

-  elastic  and  plastic  strain  increments  at  the  internal 


wal  1 


-  internal  pressure  and  wall  closures  (internal  and 
externa  1 ) 

-  radial  stress  at  the  boundary  from  both  the  elastic 
and  plastic  solutions 

-  radius  of  the  plastic  zone 

Typical  radial  and  tangential  stress  distributions 
within  the  cylinder  are  shown  in  Figure  4.6. 
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Figure  4.6*  Radial  and  Tangential  stress  distributions 

within  a  hollow  cylinder 
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4.5  Sensitivity  Analysis 

The  model  was  used  to  check  the  sensitivity  of  the 
solution  to  the  various  input  parameters  and  the  size  of  the 
iteration  grid  )  .  The  strength  parameters  (a,b)  and 
dilation  factor  (oc)  all  influence  the  solution.  The 
closures  determined  from  the  analyses  have  been  normalized 
to  the  elastic  closures.  This  normalization  procedure 
eliminates  the  influence  of  the  elastic  deformation 
parameters  on  the  solution. 

The  results  of  the  analyses  are  presented  in  Figures 
4.7  to  4.11. 

4.6  Discussion 

The  size  of  the  iteration  grid  (^\r)  only  has  a  large 
influence  on  the  closure  when  the  dilation  parameter  is 
large.  This  case  is  close  to  that  of  the  associated  flow 
rule  and,  as  has  been  discussed  before,  is  not  applicable  to 
drained  behaviour  of  sands.  It  can  thus  be  concluded  that 
the  size  of  the  iteration  grid  does  not  have  a  significant 
influence  on  the  closure. 

The  spacing  of  the  nodes  should  be  selected  with  some 
judgement  since  the  value  of  the  plastic  radius  calculated 
will  lie  between  two  values  which  differ  by  an  amount  . 
Accuracy  of  the  solution  will  depend  on  whether  the  nodes 
are  finely  or  coarsely  spaced.  Actual  variation  of  /\  r  in  a 
finely  spaced  nodal  system  will  not  have  much  influence  on 
the  closure. 
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NORMALIZED  CLOSURE 

Figure  4.7  Sensitivity  of  the  number  of  nodes  on  the 

solution 
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Figure  4.8  Sensitivity  of  the  Alpha-parameter  on  the 

solution 
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NORMALIZED  CLOSURE 

Figure  4.9  Sensitivity  of  the  a-parameter  on  the  solution 

(b=0 . 8 ) 
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NORMALIZED  CLOSURE 

Figure  4 .  10  Sensi t i vi ty  of  the  a-parameter  on  the  solution 

( b= 1.0) 
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NORMALIZED  CLOSURE 

Fiqure  4.11  sensitivity  or  the  b-parameter  on  the  solution 

( a=0 . 7 ) 
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The  dilation  parameter  (  o( )  has  a  significant  influence 
on  the  solution,  especially  when  it  has  a  value  approaching 

that  required  to  represent  the  associated  flow  rule.  The 

% 

value  of  ^  cannot  be  determined  from  laboratory  tests  and 
must  be  selected  intuitively. 

The  associated  flow  rule  is  probably  only  valid  for 
drained  behaviour  near  peak  strength,  and  behaviour  of  the 
material  after  peak  will  best  be  represented  by  a 
non-associ ated  flow  rule.  The  larger  the  plastic  strains, 
the  lower  will  the  -value  be.  The  plastic  strains  in  the 
hollow  cylinder  test  are  large,  and  a  value  of  ^  =1  (i.e.  no 
volume  change)  will  probably  be  the  best  value  to  use  in  the 
analysis.  This  value  of  o4  was  used  in  the  analysis  of  the 
Ottawa  sand  sample  (see  Figure  4.5)  when  comparing  predicted 
closures  with  the  measured  ones. 

The  strength  parameters  (a,b)  have  two  influences  on 
the  results.  The  higher  the  value  of  either  of  these 
parameters,  the  higher  is  the  yield  point  (the  stress  at 
which  yielding  first  occurs).  The  yield  point  does  not  seem 
to  be  significantly  affected  by  a  small  change  in  the 
a-value  when  b  is  less  than  unity. 

The  closures  are  also  affected  by  the  strength 
parameters.  The  higher  the  value  of  a  or  b,  the  smaller  are 
the  closures.  This  is  to  be  expected  since  the  comparisons 
are  made  at  equal  internal  pressures.  In  a  stronger 
material,  less  plastic  flow  has  occured  at  the  same  pressure 


ratio. 
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In  Figure  4.8,  the  results  shown  in  the  centre  diagram 
are  completely  different  to  the  other  two  diagrams.  These 
results  have  been  rechecked,  and  no  error  was  found.  A 
slight  variation  in  the  a  or  b  values  resulted  in  a  set  of 
curves  similar  to  the  other  two.  It  must  be  concluded  that 
the  analytical  procedure  is  unstable  for  this  particular  set 
of  data.  The  effects  of  grid  size  and  size  of  the  stress 
increments  may  influence  this  instability. 

The  variation  of  E  in  the  elastic  region  of  the 
cylinder  must  be  correctly  specified.  A  constant  E  value 
will  result  in  a  large  under-predict ion  of  the  deformations. 
The  value  of  E  is  a  function  of  the  degree  of  confinement 
within  the  elastic  region  as  well  as  the  stress  level.  The 
variation  of  E  within  the  elastic  region  can  be  modelled  by 
a  bi-linear  variation  of  the  modulus  with  the  average  radial 
stress,  as  determined  from  the  hollow  cylinder  test. 

The  results  of  the  sensitivity  analysis  have  shown  the 
importance  of  careful  selection  of  the  input  parameters. 

The  values  of  the  strength  parameters,  dilation  parameter 
and  the  nodal  spacing  can  have  a  large  influence  on  the 
predicted  deformations.  While  the  magnitude  of  Young's 
modulus  will  affect  the  closures,  a  more  detrimental  effect 
will  result  if  the  variation  of  modulus  with  the  degree  of 
confinement  is  not  considered. 

The  model  has  been  developed  for  the  purpose  of 
predicting  the  deformations  that  occur  in  the  hollow 
cylinder  test.  Although  this  has  been  the  main  objective, 


. 
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the  model  can  also  be  used  to  analyse  the  deformations  that 
occur  in  deep  tunnels  or  around  shafts  (i.e.  the  stress 
field  must  be  isotropic). 

The  plane-strain  solution  is  applicable  for  tunnel 
analysis,  but  the  gravitational  effects  must  still  be 
considered.  A  shaft  analysis  is  a  plane-stress  problem,  but 
it  has  been  suggested  by  McCreath ( 1 98 1 )  that  the  problem  can 
be  analysed  with  sufficient  accuracy  by  the  plane-strain 
solution.  Gravitational  effects  do  not  influence  the 
solution  for  the  deformations  that  occur  around  a  shaft. 

The  model  has  been  shown  to  predict  the  closures  in  a 
hollow  cylinder  test  with  sufficient  accuracy  (see  also 
Chapter  6).  The  finite  boundaries  of  the  cylinder  make  the 
use  of  determining  the  elastic  modulus  from  the  average 
radial  stress  in  the  elastic  region  (Assumption  #1)  viable. 
In  the  case  of  a  shaft  or  tunnel,  the  external  boundary  will 
not  be  subjected  to  any  significant  movements  and  hence  the 
use  of  Assumption  #3  to  determine  the  elastic  modulus  will 
probably  yield  a  better  solution.  It  was  shown  in  Figure 
4.5  that  this  assumption  actually  gave  a  better  prediction 
of  the  internal  wall  movement,  but  did  not  give  a  good 
prediction  of  the  external  wall  movement. 

The  model  may  be  modified  to  suit  any  particular 
requirement.  The  determination  of  the  elastic  modulus  may 
be  found  by  any  method  with  relatively  few  changes  required 
in  the  program.  The  failure  criterion  may  also  be  specified 
differently  if  required.  The  criterion  used  can  incorporate 
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a  strength  loss  (strain-weakening)  with  relative  ease.  This 
feature  has  not  been  included  in  the  model  since  the 
deformations  of  the  cylinder  can  be  predicted  reasonably 
well  without  it.  The  parameters  that  describe  the  strength 
loss  cannot  be  measured  in  the  test  and  this  results  in 
another  variable  that  has  to  be  experimented  with  during  the 
analysi s . 

A  non-linear  elast ic-perfect ly  plastic  model  has  thus 
been  developed  which  can  predict  the  deformations  that  occur 
in  the  hollow  cylinder  test  with  reasonable  accuracy.  This 
model  is  an  analytical  tool  that  will  assist  in  the  data 
interpretation  of  these  tests.  It  also  has  the  possibility 
of  being  used  in  the  analysis  of  the  deformations  that  occur 
around  shafts  and  tunnels  at  fairly  deep  locations. 


CHAPTER  5 


EXPERIMENTAL  RESULTS 


5.1  Introduction 

The  influence  of  stress  path  on  the  deformation 
behavior  of  oil  sands  is  examined.  A  laboratory  program  was 
conducted  in  which  oil  sands  were  tested  under  three 
different  stress  paths.  The  triaxial  testing  of  the  samples 
was  conducted  in  a  hollow  cylinder  triaxial  device  and  in  a 
high  pressure  triaxial  cell. 

An  oedometer  test  was  conducted  to  determine  the 
compressibility  of  the  soil  skeleton,  and  to  compare  this 
result  with  those  determined  from  the  triaxial  tests.  Index 
properties  of  the  samples  were  determined. 

The  material  tested  was  an  oil  rich,  dense  sand  with 
cross-bedding.  Sample  collection  and  preparation  techniques 
have  been  presented  in  Chapter  2.  All  samples  were  taken 
perpendi cu 1 ar  to  the  cross-bedding. 


5.2  Laboratory  Procedures 
5.2.1  Index  Tests 

Index  tests  were  conducted  to  determine  grain  size 
distribution,  density,  porosity  and  fluid  content  (bitumen 
and  water).  Grain  size  distributions  were  conducted  on  both 
the  original  samples  and  the  samples  after  completion  of  the 
tests.  The  samples  were  washed  on  a  number  200  U.S.  sieve 
(0.074  mm)  and  the  analysis  conducted  according  to  ASTM 
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standards.  The  material  passing  the  #200  sieve  has  been 
designated  as  the  percentage  fines  of  the  sample. 

The  densities  were  determined  from  the  prepared 
samples,  whose  dimensions  and  weight  were  recorded  prior  to 
testing.  Fluid  contents  were  obtained  by  conducting 
extraction  tests  on  the  samples  after  the  test.  Original 
fluid  contents  were  estimated  from  representat i ve  material 
taken  from  the  block  samples. 

5.2.2  Oedometer  Test 

The  testing  apparatus  used  was  designed  and  built  at 
the  University  of  Alberta.  The  frame  had  the  capability  of 
applying  normal  stresses  of  up  to  16  MPa  to  the  sample  by 
means  of  dead  weights  loaded  onto  a  lever  arm.  The  applied 
loads  were  measured  on  a  44  kN  load  cell  and  the  vertical 
displacements  on  a  linear  variable  differential  transducer. 
The  compressibility  of  the  soil  skeleton  was  determined  from 
the  slope  of  the  volumetric  strain-effective  stress  curve. 

5.2.3  Triaxial  Tests 

The  test  procedures  and  sample  preparation  for  the 
hollow  cylinders  and  triaxial  samples  have  been  deer i bed  in 
detai 1  in  Chapter  2 . 
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5.3  Laboratory  Results 

5.3.1  Index  Tests 

5. 3. 1.1  Grain  Size  Distribution 

The  grain  size  distribution  for  all  samples  are  given 
in  Appendix  E  for  both  the  original  samples  and  for  those 
after  completion  of  the  tests.  A  typical  grain  size  curve 
of  the  oil  sands  tested  is  given  in  Figure  5.1.  It  can  be 
seen  that  the  material  is  a  uniform,  fine-grained,  micaceous 
sand  with  a  small  percentage  of  fines. 

An  increase  in  the  amount  of  fines  after  completion  of 
the  oedometer  and  triaxial  tests  indicates  that  grain 
crushing  occurs  under  high  pressures.  Table  5.1  shows  the 
results  of  the  grain  size  analyses. 

5. 3. 1.2  Density  and  Fluid  Contents 

The  results  of  these  tests  are  presented  in  Table  5.2. 
The  bulk  density  of  the  samples  range  from  1.93-1.99  g/cc. 
This  density  is  higher  than  would  be  expected  for  a  dense 
sand  and  is  indicative  of  the  interlocking  nature  of  the 
sand . 

The  initial  fluid  content  (bitumen  +  water)  of  1 7% - 19% 
is  equivalent  to  a  degree  of  saturation  of  approximately 
80%.  This  value  is  low  and  due,  to  some  degree,  to  the 
drying  of  the  blocks.  However,  the  low  saturation  is  also 
due  to  expansion  of  the  sample.  The  initial  porosity  of  the 
blocks  were  approximately  35%,  which  is  higher  than  that 
expected  in  situ. 
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Table  5.1  Results  of  the  grain  size  analyses 


Sample  # 

Press . 

(MPa) 

V.  Fines 

Passing  No.  200  sieve 

Before 

After 

Oedometer 

16.0 

20.4 

23.4 

A-  1 

7.3 

2.6 

6.4 

CM 

i 

<t 

10.0 

2.6 

4  .  9 

B-1 

7.4 

6.7 

8.  1 

CM 

CD 

10.1 

6.7 

12.  1 

C-1 

7.0 

6.7 

10.  1 

02 

4.3 

6.7 

5.4 

Table  5.2  Results  of  the  Index  tests 
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Figure  5.1  Typical  grain  size  distribution  for  oil  sand 

samples 
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5. 3. 1.3  Oedometer  Results 

The  compressibility  of  oil  sands  is  stress  dependent, 
as  can  be  seen  from  Figure  5.2.  The  compressibility  remains 
relatively  constant  after  one  load-unload  cycle,  and  it  is 
not  necessary  to  conduct  cyclic  compressibility  tests  on  the 
triaxial  samples  prior  to  shearing.  The  results  of  the  test 
are  shown  in  Table  5.3(a).  The  compressibility  of  the  oil 
sands  tested  is  higher  than  that  reported  by  Barnes  (1980) 
and  similar  to  those  reported  by  Dusseault  (1980).  These 
higher  compressibilities  indicate  that  some  sample 
disturbance  has  occurred.  The  reported  compressibilities 
are  given  in  Table  5.3(b)  and  (c). 

5.3.2  Triaxial  Tests 

5.3.2. 1  Strength  Results 

The  results  of  the  compressibility  test,  prior  to 
shearing,  of  the  samples  are  given  in  Appendix  H,  and  the 
data  shown  in  Table  5.3.  The  stress-strain  curves  obtained 
from  the  three  series  of  tests  are  shown  in  Figures  5.3  to 
5.5  along  with  the  volumetric  strain-effective  stress 
curves.  It  should  be  noted  that  the  stress-strai n  curves 
for  Series  A  (hollow  cylinder  tests)  are  based  on  the 
assumption  of  a  linear  radial  stress  distribution  within  the 

sample . 

The  axial  stress  vs.  radi a  1  +tangent i a  1  stress  results 
for  Series  A  are  shown  in  Figures  5.6(a)  and  5.6(b).  These 
curves  are  used  to  determine  Poisson  s  ratio  for  the 
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Table  5.3  Compressibility  data 


(a)  Test  results 


Sample  # 

Compressibility,  Cs 
(kPa-1  )  *  ]0~6 

4  MPa 

7  MPa 

10  MPa 

Oedometer 

1st  Loading 

3.46 

2.60 

2.60 

Unloading 

0.90 

0.61 

Reload  i  ng 

0.65 

0.65 

2nd  Unloading 

0.22 

0.22 

A-1 

1.46 

1.06 

A-2 

1.61 

0.96 

0.55 

B-1 

3.35 

2.56 

B-2 

4.12 

2.82 

2.07 

C  - 1 

3.60 

2.64 

C— 2 

3.46 

(b)  After  Barnes  (1980) 


Material 

Modulus  of  Compressibility,  av  (kPa*'*10~6) 

Load 

Unload 

Reload 

Unload 

Fine-  grained 

Me  Murray  Formation 
(oil  -  free ) 

0.80 

0,33 

0.35 

0.28 

Note:  Cs  =  =  ~fV  (for  n=32V,) 


(c)  After  Dusseault  (1980) 


Material 

Compressibility 
(kPa"' )  »  10'* 

Highly  disturbed  samples 

High  quality  core 

Probable  in  situ  value 

10'3  to  10"5 

5  x  10"6 

10"6  to  10'7 

' 

VOL.  STRAIN  m 
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Figure  5.2  Oedometer  compressibility  test 
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Figure  5.3  Stress-strain  curves  for  Series  A 
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Figure  5.4  Stress-strain  curves  for  Series  B 
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Figure  5.5  Stress-strain  curves  for  Series  C 
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Figure  5.6  Axial  stress  vs.  radi a  1 +  tangent i a  1  stress 

(Series  A) 


88 


p 1 ane-strai n  condition  of  these  tests.  The  effective  stress 
paths  for  all  tests  are  shown  in  Figure  5.7,  and  the  results 
of  the  tests  are  summarized  in  Table  5.4.  The  effective 
friction  angle  (  )  is  defined  as  the  secant  angle  of 

friction: 


The  angle  of  friction  is  calculated  from  the  stresses 
at  failure.  For  the  triaxial  tests  (Series  B  and  C), 
failure  is  defined  as  the  maximum  deviator  stress.  The 
stresses  at  failure  for  the  hollow  cylinder  tests  (Series  A) 
are  defined  as  the  stress  state  just  prior  to  collapse  of 
the  sample. 

The  failure  envelope  for  Series  B  and  C  tests  is  shown 
in  Figure  5.8.  It  can  be  seen  that  the  envelope  is 
approximately  linear  in  the  stress  range  considered,  but 
does  not  pass  through  the  origin.  This  indicates  that 
curvature  of  the  envelope  would  occur  in  the  lower  stress 
range,  since  the  oil  sands  have  no  shear  strength  at  zero 
effective  stress  ( Dusseaul t , 1 978 ) . 

The  results  of  Series  B  and  C  tests  have  been  compared 
with  similar  tests  conducted  by  Dusseault  (1978,1981)  on 
both  Ottawa  sand  and  oil  sand.  The  comparison  is  shown  in 
Figure  5.9.  This  comparison  is  useful  in  analysing  the 
results.  Dusseault' s  results  are  obtained  from  both  lean 
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Table  5.4  Summary  of  test  results 
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Figure  5.7  Effective  stress  paths  for  all  tests 
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Figure  5.8  Strength  envelope  for  Series  B  and  C  tests 
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Figure  5.9  Comparison  of  test  results  with  published  data 
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and  rich  oil  sand  samples  of  different  diameters.  The 
results  indicate  that  the  strength  of  the  oil  sand  as 
measured  in  this  laboratory  program  are  much  lower  than  the 
reported  strengths. 

Discussion  of  these  results  is  given  in  a  later 
section . 


5. 3. 2. 2  Deformation  Results 

The  influence  of  stress  path  on  the  elastic  deformation 
moduli  can  be  examined  in  three  ways.  The  variation  of  the 
Young's  modulus  with  consolidation  pressure  for  each  series 
of  tests  can  be  compared  (Figure  5.10(a)),  the  volumetric 
strain  at  failure  can  be  compared  at  different  consolidation 
pressures  for  the  three  stress  paths  (Figure  5.10(b)),  or 
the  Poisson's  ratio  can  be  compared  (Figure  5.10(c)).  When 
comparison  of  the  Young's  moduli  are  made,  it  is  convenient 
to  define  the  modulus  as  the  slope  of  the  stress-strain 
curve  between  two  particular  stress  levels.  Two  moduli  are 
defined  in  this  manner  as  given  below: 

E-  - 
1 


Eeff 


^50  ~  °25 

^50  "  ^25 


where  C^Si^O  are  stresses  at  25V.  and  50V.  ot  the 

maximum  deviator  stress. 

625 ,  £50  are  strains  corresponding  to  O25  and  O50 
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Figure  5 . 1 0  Inf luence  of  the  stress  path  on  the  deformation 

behaviour 
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While  this  approach  may  not  give  the  actual  modulus  at 
any  particular  stress  level,  it  makes  comparisons  simpler. 
The  deformation  moduli  are  summarized  in  Table  5.5. 

Comparison  of  these  plots  show  that  the  deformation 
moduli  are  dependent  on  the  stress  path  followed  in  the 
test.  These  results  are  discussed  in  more  detail  later. 


5.4  Failure  Modes  of  Triaxial  Samples 

5.4.1  Hollow  Cylinder  Tests  (Series  A) 

The  samples  were  carefully  disassembled  after  testing. 
The  hollow  cylinders  still  had  an  approximately  circular 
cross-section  and  were  self-supporting.  They  were  then 
frozen  in  a  cold  room,  sectioned,  and  photographed. 

Plate  5.1  shows  the  condition  of  Samples  #A-1  and  #A-2 
immediately  after  testing.  The  failure  surface  is  clearly 
visible  for  Test  #A-2  as  the  membrane  is  distorted  in  the 
vicinity  of  the  shear  plane. 

The  ha  1 f -sect ions  of  the  frozen  samples  are  shown  in 
Plates  5.2  and  5.3  for  Tests  #A-1  and  #A-2  respectively. 

The  failure  surface  for  Test  #A-1  is  not  distinctive  but  a 

* 

trace  of  the  surface  can  be  seen  on  the  inside  and  outside 
walls  of  the  cylinder.  This  failure  surface  is  only  found 
in  one  location  i.e.  there  is  no  shearing  across  the  entire 
sample.  The  failure  mode  is  typical  for  the  hollow  cylinder 
test  when  sheared  under  such  conditions.  Sample  #A-2  failed 
in  an  identical  manner,  but  the  failure  surface  is  much  more 


di st i net . 


. 


' 
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Table  5.5  Summary  of  the  deformation  moduli 


Sample  # 

Ini t i al 

Tangent 

Modulus 

( GPa ) 

Eff . 

T  angent 

Modu 1  us 

(GPa) 

Poisson' s 

Ratio 

A-  1 

3.50 

3.00 

0.20-0.21 

A-2 

3.00 

2.00 

0.20-0.21 

B-1 

0.982 

0.872 

0.24-0.26 

B-2 

0.776 

0.593 

0.13-0.17 

C-1 

0.848 

0.714 

0.37-0.42 

C-2 

0.496 

0.419 

0.37-0.42 
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Plate  5.1  Condition  of  samples  A-1  and  A-2  immediately  after 

the  tests 
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Plate  5.2  Ha  1 f -sect ions  of  sample  A-1  after  the  test 
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Plate  5.3  Half-sections  of  sample  A-2  after  the  test 
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The  actual  shape  of  the  failure  surface  through  the 
wall  of  the  sample  cannot  be  examined  due  to  the  bitumen  in 
the  material.  However,  from  the  trace  of  the  surfaces  on 
both  inside  and  outside  walls,  it  is  likely  that  the  failure 
surface  is  spiral  in  plan  view.  Failure  occurs  along  the 
most  kinematically  possible  surface. 

The  trace  of  this  surface  is  inclined  at  approximately 
70  degrees  to  the  horizontal  in  the  middle  third  of  the 
sample  height.  This  surface  then  changes  direction  at  the 
top  of  the  sample  due  to  the  end  effect  of  the  loading  cap. 
The  end  effect  of  the  bottom  of  the  sample  is  not  as 
pronounced . 

Examination  of  the  half-section  of  the  cylinder  showed 
the  bore  to  have  deformed  in  a  uniform  manner.  No  deviation 
from  uniform  movement  could  be  seen,  even  near  the  ends  of 
the  sample.  This  uniform  deformation  supports  the  use  of  an 
axi -symmetr ic  analysis  of  the  sample,  as  well  as  the  use  of 
volume  change  measurements  to  determine  the  strains  that 
occur  in  the  sample  (which  are  assumed  to  be  uniform  along 
the  height  of  the  cylinder). 

5.4.2  Triaxial  Tests  (Series  B  and  C) 

At  the  completion  of  the  test,  the  samples  were  removed 
from  the  membrane  and  examined.  No  distinct  shear  planes 
were  visible  and  the  overall  deformation  of  the  sample  was 
uniform.  No  bulging  of  the  sample  occurred. 

The  end  restraint  of  the  load  cap  resulted  in  a  minor 
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distortion  at  the  top  of  the  sample. 


5.5  Discussion  of  Results 

The  strength  of  the  hollow  cylinder  tests  depend  on  the 
analysis  used,  and  hence  have  not  been  included  in  Figure 
5.9.  A  detailed  analysis  of  this  series  of  tests  is 
presented  in  Chapter  6. 

The  strength  envelopes  shown  in  Figure  5.9  reflect  the 
influence  of  bitumen  content  on  the  strength  of  oil  sand. 

The  data  from  lean  oil  sands  tested  by  Dusseault  (1978) 
provide  an  upper  bound  to  the  shear  strength  while  the 
strengths  determined  from  this  testing  program  provide  a 
lower  bound.  The  bitumen  per  se  does  not  contribute  to  the 
strength  of  oil  sands  ( Dusseaul t , 1 978 )  but  has  a  dominant 
role  in  the  amount  of  sample  disturbance  that  will  occur. 

For  the  rich  oil  sand,  the  exsolution  of  gas  that 
occurs  when  unloading  the  material  (during  sampling)  results 
in  an  expansion  of  the  material  and  a  disruption  of  the 
fabric.  This  results  in  a  significant  drop  in  the  shear 
strength . 

The  original  porosity  of  the  samples  (35%)  is  higher 
than  the  in  situ  porosity  but  still  reflects  a  porosity 
lower  than  for  most  dense  sands.  A  comparison  of  the 
porosities  immediately  after  preparation  of  the  samples 
shows  an  increase.  This  increase  is  larger  for  the  triaxial 
samples  ( n=37 . 6%-38 . 5% )  than  for  the  hollow  cylinder  ones 
(n=36%).  These  porosity  changes  are  summarized  in  Table 
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5.6. 

The  increase  in  porosity  after  sample  preparation  is 
due  to  sample  disturbance,  and  is  reflected  more  in  a 
smaller  sample  than  in  a  larger  one.  The  dissipation  of  the 
heat  generated  when  trimming  a  sample  will  depend  on  the 
volume  of  the  sample.  The  small  tri axial  samples  (37.5mm) 
have  a  smaller  volume  and  therefore  heat  up  more.  This 
results  in  more  expansion  as  is  evident  from  the  porosity 
measurements.  During  trimming  of  the  hollow  cylinder 
samples,  the  bore  was  Kept  cold  by  using  dry  ice  in  the 
shaft.  This  helped  in  dissipating  the  heat  generated  by 
tr immi ng . 

The  tests  on  rich  oil  sands  conducted  by  Dusseault 
(1978)  were  carried  out  at  lower  stress  levels  on  larger 
diameter  samples  (75  mm).  These  results  show  a  large 
curvature  of  the  envelope  and  tend  to  correspond  with  the 
test  results  of  this  program.  One  data  point  from  those 
tests  coincides  with  the  envelope  obtained  for  lean  oil 
sands.  This  indicates  that  the  oil  sand  has  a  higher 
strength  than  the  majority  of  the  tests  would  indicate.  The 
reduction  in  ^strength  is  due  again  to  disturbance  of  the 
material  during  sampling  and  possibly  during  preparation  of 
the  samples. 

The  rich  oil  sand  samples  tested  yield  lower  strengths 
than  obtained  for  Ottawa  sand.  The  oil  sand  is  a 
fine-grained,  sub-angular  sand  and  should  thus  have  a  higher 
strength  than  the  coarser-grai ned ,  rounded  Ottawa  sand. 
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Table  5.6  Influence  of  sample  disturbance  on  porosity  and 
densi ty 


MATERIAL 

it 

wbit 

C/.) 

Before 

Preparation 

After 

Preparation 

& 

(g/cc ) 

n 

(%) 

(g/cc  ) 

n 

(%) 

Block  used  for 
Oedometer  sample 

2.0 

(1.7) 

2.01 

35.2 

2.00 

35.5 

Block  used  for 
Series  B  &  C 

16.6 
(14.1  ) 

2.01 

35.5 

1.93 

-1.96 

34.6 

-38,5 

Block  used  for 
Series  A 

16.2 

(13.7) 

2.01 

35.2 

1.98 

-1.99 

36.0 

-36,8 

tt  Figures  in  brackets  based  on  sample  volume 


rv 
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However,  the  initial  porosity  of  the  Ottawa  sand  is 
approximately  34%,  while  that  of  the  oil  sand,  after 
consolidation,  is  only  35%  to  36%.  This  difference  in 
porosity  will  have  a  significant  effect  on  the  frictional 
resistance  of  the  material.  Variations  in  grain  size  and 
shape  will  also  have  some  influence  on  the  shear  strength, 
but  it  is  the  difference  in  porosities  that  will  have  the 
largest  effect. 

The  results  show  that  once  an  oil  sand  sample  is 
disturbed,  its  original  fabric  is  destroyed  and  will  not  be 
recovered  by  subsequent  reconsolidation.  Consolidation  of  a 
sub-angular  sand,  even  to  pressures  of  12  MPa,  will  not 
result  in  a  porosity  that  can  be  obtained  by  vibrating  a 
rounded  sand.  It  is  thus  extremely  important  to  minimize 
sample  disturbance  of  oil  sands. 

The  drained  strength  of  rich  oil  sand  has  not  been 
determined  with  any  degree  of  confidence.  A  lower  bound 
value  from  the  test  results  may  be  used,  but  this  will 
result  in  conservative  design  analysis.  It  has  been  shown 
that  the  careful  and  elaborate  procedures  used  to  prepare 
samples  of  low  gas  saturation  has  not  been  sufficient  to 
prevent  sample  disturbance.  This  will  pose  a  more  serious 
problem  when  trying  to  prepare  gas  saturated  samples 
obtained  at  depth. 

The  laboratory  preparation  techniques  will  have  to  be 
improved.  Coring  of  hollow  cylinder  samples  can  be  improved 
by  using  a  drilling  fluid  that  has  been  cooled,  but  such  a 
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fluid  must  not  penetrate  the  walls  of  the  sample.  Trimming 
techniques  provide  the  most  problems,  and  a  system  for 
reducing  the  heat  generation  needs  to  be  examined. 

The  strains  at  failure  of  the  tests  also  indicate  the 
effect  of  sample  disturbance.  In  Series  A,  the  strain  to 
failure  ( 1 . 3% - 2% )  is  typical  for  a  brittle  material  while 
those  for  Series  B  and  C  ( 4% - 7 % )  are  more  typical  of  a  dense 
sand  at  high  pressures  (ductile).  Same  of  this  difference 
may  be  due  to  the  plane-strain  condition  for  the  hollow 
cylinder  tests. 

The  influence  of  stress  path  on  the  deformation  moduli 
have  been  shown  in  Figures  5.10(a)  to  5.10(c).  The  data  are 
not  sufficient  to  make  direct  conclusions,  but  the  trends 
shown  do  indicate  that  deformation  moduli  are  path 
dependent.  Young's  modulus  for  all  three  series  of  tests 
show  different  behavior.  The  hollow  cylinder  tests  yield 
higher  moduli  but  this  could  be  due  to  the  relative  amounts 
of  sample  disturbance.  However,  the  trend  of  Series  B  and  C 
are  different  even  though  the  sample  preparation  and  initial 
porosities  are  the  same. 

The  Poisson's  ratio  is  dependent  on  the  failure  mode  of 
the  material.  Series  A  and  C  show  no  dependency  on  the 
consolidation  pressure.  In  Series  B,  the  modulus  decreases 
with  consolidation  pressure.  This  behaviour  is  a  result  of 
the  different  failure  modes  that  occur  in  this  series  of 
tests  (dilation  and  contraction).  The  influence  of  stress 
path  can  be  best  examined  from  Figure  5.10(b).  The  stress 
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path  followed  in  Series  A  and  C  (see  Chapter  2)  are  similar, 
and  their  volume  change  behaviour  is  also  similar.  The 
volume  change  behaviour  of  Series  B  is  completely  different. 
At  a  higher  consolidation  pressure  (lOMPa),  dilatancy  ts 
suppressed  while  at  a  pressure  of  7  MPa  the  volumetric 
strain  is  more  negative  (dilatant)  than  for  the  other  two 
ser i es . 

The  results  shown  in  Figure  5.10  show  the  dependency  of 
the  deformation  behaviour  on  stress  path.  However,  the 
material  has  been  disturbed  and  these  trends  may  be 
completly  different  for  the  intact  material.  The  difference 
between  Series  B  and  C  should  still  be  the  same.  For  an 
undisturbed  sample,  suppression  of  dilatancy  will  probably 
occur  at  a  higher  confining  pressure  while  the  J1 -constant 
test  at  this  consolidation  pressure  still  requires  a 
reduction  of  cell  pressure  during  the  test,  and  hence  will 
not  suppress  dilation. 


5.6  Conclusions 

The  influence  of  stress  path  on  the  deformation 
behaviour  of  oil  sands  must  be  considered  in  any  analysis. 
The  stress  path  associated  with  deformation  around  tunnels 
or  shafts  can  be  approximated  in  a  hollow  cylinder  triaxial 
test  or  in  a  Jl-constant  triaxial  test.  A  triaxial 
compression  test  is  not  relevant  to  such  analyses  and  will 
provide  incorrect  deformation  properties. 
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The  results  of  the  laboratory  tests  show  the  influence 
of  the  sample  preparation  techniques  and  sample  collection 
on  the  disturbance  of  rich  oil  sands.  These  procedures  need 
to  be  perfected  before  reliable  strength  and  deformation 
properties  can  be  determined. 
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CHAPTER  6 


PREDICTION  OF  THE  HOLLOW  CYLINDER  RESULTS  FROM  THE  NUMERICAL 

MODEL 


6.1  Introduction 

The  deformation  behaviour  of  the  oil  sands  as 
determined  from  the  hollow  cylinder  tests  was  given  in  the 
previous  chapter.  Those  results  were  based  on  the 
assumption  of  a  linear  radial  stress  distribution  within  the 
•cylinder.  In  this  Chapter,  the  deformation  behaviour  is 
examined  by  comparing  stress-strain  curves  for  the  tests 
based  on  different  assumptions. 

The  closures  measured  in  the  tests  are  compared  with 
those  predicted  by  the  numerical  model  developed  in  Chapter 
4.  The  model  is  also  used  to  determine  the  strength  of  the 
cylinders  and  the  values  are  compared  with  the  actual  test 
data.  The  influence  of  the  assumptions  made  (regarding 
stress  distribution)  when  determining  the  cylinder  strength 
from  the  test  data  is  examined. 


6.2  Analytical  Procedure 

The  procedure  adopted  in  the  analysis  is  shown  in  the 
flow  chart  in  Figure  6.1.  The  test  data  is  used  to  estimate 
the  strength  parameters  and  to  determine  the  relationship 
between  Young's  modulus  and  average  radial  stress.  Values 
for  the  di lation  parameter  (  cK)  and  Poisson  s  ratio  (  V  )  are 

selected . 
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Figure  6.1  Flow  chart  showing  the  analytical  procedure 
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These  values  are  used  as  the  initial  input  parameters. 
The  strength  parameters  are  adjusted  until  the  cylinder 
collapses  at  the  same  pressure  ratio  as  occurs  in  the  actual 
test.  The  predicted  closures  are  compared  with  the  actual 
closures  and  the  input  parameters  can  be  adjusted  until  the 
best  comparison  is  achieved. 

The  strength  parameters  (a,  b)  required  to  achieve 
equivalent  collapse  of  the  cylinder  are  compared  with  those 
determined  from  the  test  data.  These  parameters  are  used  to 
assess  the  validity  of  the  different  assumptions  made 
regarding  the  stress  distribution  within  the  cylinder. 


6.3  Analysis  of  the  Hollow  Cylinder  test 

Two  assumptions  were  made  when  determining  the  average 
radial  stress  in  the  cylinder.  The  first  was  that  the 
stress  varied  linearly  throughout  the  cylinder,  while  the 
second  assumption  was  based  on  plasticity  theory  to 
determine  the  stress  distribution.  The  theory  was 
simplified  by  assuming  that  the  value  of  (see  Chapter  3) 
after  yielding  occurs  is  the  same  as  the  value  of  (JQ  with 
no  yielding. 

The  two  hollow  cylinder  tests  conducted  on  oil  sand 
samples  as  well  as  one  conducted  on  an  Ottawa  sand  sample 
are  examined  in  the  following  sections. 


■ 
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6.3.1  Sample  #A- 1 

The  stress-strain  curves  for  both  assumptions  are  shown 
in  Figure  6.2,  and  the  measured  closures  are  given  in  Figure 
6.3  for  both  the  internal  and  external  walls.  The  E-(CTr)ave 
relationship  determined  from  the  stress-strain  curve  is 
shown  in  Figure  6.4. 

The  failure  criterion  was  based  on  a  linear 
Mohr-Coulomb  envelope  (i.e  b= 1.0)  since  the  secant  angle  of 
friction  is  determined  from  actual  test  data  on  this  basis. 

The  comparison  of  the  predicted  closures  with  the 
measured  ones  is  shown  in  Figure  6.5.  The  predicted 
closures  are  based  on  an  assumed  oC  -value  of  1.0  (no  plastic 
volume  change  occurring)  and  a  Poisson's  ratio  of  0.25. 

The  strength  parameter  (a)  determined  from  the  analysis 
corresponds  to  a  secant  angle  of  friction  which  passes 
through  the  top  of  the  Mohr  circle.  This  value  is  then 
converted  to  a  secant  friction  angle  which  is  tangent  to  the 
Mohr-circle.  This  allows  for  comparisons  to  be  made  with 
the  measured  values  and  the  theoretical  value  as  determined 
from  Equation  3.14  (after  Wu  et  a/., 1963).  The  comparison 
of  the  strength  parameter  is  given  in  Table  6.1  along  with 
the  other  two  tests. 

6.3.2  Sample  #A-2 

The  stress-strain  curves  and  measured  closures  from  the 
test  are  given  in  Figures  6.6  and  6.7.  The  variation  of  E 
with  average  radial  stress  is  shown  in  Figure  6.8.  The 
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Table  6.1  Strength  parameters  determined  from  the  different 

assumptions 
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Figure  6.2  Stress-strain  curves  :  Test  #A-1 
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Figure  6.3  Measured  wall  movements  :  Test  #A-1 
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Figure  6.5  Comparison  of  predicted  closure  with  actual 

closure  :  Test  #A-1 
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Figure  6.6  Stress-strain  curves  :  Test  #A-2 


PKKPA)  XI  O'* 


118 


Figure  6.7  Measured  wall  movements  :  Test  #A-2 
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Figure  6.8  Variation  of  E  with  average  radial  stress  :  Test 
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predicted  closures  are  compared  with  the  actual  closures  in 
Figure  6.9,  and  the  strength  parameters  are  shown  in  Table 
6.1.  The  Young's  moduli  used  have  also  been  altered  to 
include  the  variation  of  E  in  the  very  low  stress  range  (see 
Figure  6.6).  These  results  are  shown  in  Figure  6.10. 

6.3.3  Ottawa  Sand  Sample 

This  test  has  been  included  to  compare  with  the 
closures  that  occur  in  oil  sand  samples.  The  porosity  of 
the  sample  was  33%  and  the  bulk  density  1.77  g/cc.  The 
grain  size  distribution  for  the  sand  is  shown  in  Figure 
6.11. 

The  stress-strai n  curves  are  shown  in  Figure  6.12  and 
the  closure  curves  in  Figure  6.13.  The  E-  (G*r)ave 
relationship  is  shown  in  Figure  6.14,  and  the  predicted 
closures  are  compared  with  the  actual  ones  in  Figure  6.15. 
The  strength  parameters  are  shown  in  Table  6.1. 


6.4  Discussion  of  Results 

6.4.1  Shear  strength  of  the  cylinders 

The  initial  elastic  portion  of  the  stress-strai n  curves 
for  all  three  tests  are  not  significantly  influenced  by  the 
assumption  made  for  the  stress  distribution.  However,  once 
yielding  has  occured,  and  especially  when  the  plastic  radius 
is  large,  the  stress-strain  behaviour  of  the  cylinder  is 
strongly  influenced  by  the  assumption  made. 
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Figure  6.9  Comparison  of  predicted  closure  and  actual 

closure  :  Test  #A-2 
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Figure  6.10  Comparison  of  predicted  closure  and  actual 
closure  :  Test  #A-2  (adjusted  moduli) 
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Figure  6.11  Grain  size  distribution  for  Ottawa  sand 
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Figure  6.12  Stress-strain  curves  :  Ottawa  sand  sample 
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Figure  6.13  Measured  wall  movements  :  Ottawa  sand  sample 
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Figure  6.14  Variation  of  E  with  average  radial  stress  : 

Ottawa  sand  sample 
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Figure  6.15  Comparison  of  predicted  closure  and  actual 

closure  :  Ottawa  sand  sample 
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The  E- (G*r)ave  relationship  used  was  determined  from  the 
stress-strai n  curve  based  on  the  linear  radial  stress 
assumption.  As  has  been  shown,  this  curve  is  similar  to  the 
the  second  one  in  the  elastic  range,  and  it  is  this  portion 
of  the  curve  that  is  used  to  determine  the  E-  ( CTr )ave 
relationship. 

The  use  of  this  simplified  stress-s trai n  curve  is  then 
valid  when  analysing  elastic  deformation  behaviour.  The 
maximum  deviator  stress  is  underestimated  and  hence  lower 
strength  properties  are  determined. 

The  second  assumption  made  may  also  be  used  in 
determining  elastic  deformation  behaviour.  However,  the 
assumption  largely  overpredicts  the  maximum  deviator  stress 
and  hence  the  strength  parameters. 

The  strength  parameters  as  determined  by  Wu  et 
al . (1963)  are  identical  to  those  determined  from  the 
numerical  model.  This  gives  confidence  in  the  analytical 
capabilities  of  the  model. 

These  friction  angles  are  larger  than  those  determined 
from  the  37.5  mm  diameter  triaxial  samples  by  up  to  18 
degrees.  Plane-strain  tests  usually  have  a  friction  angle 
of  about  4  degrees  higher  than  conventional  triaxial  tests. 
This  large  difference  is  due  to  the  amount  of  sample 
disturbance  that  has  occurred  prior  to  testing,  as  has  been 
discussed  in  a  previous  section. 

These  results  indicate  that  the  hollow  cylinder 
triaxial  samples  are  still  relatively  undisturbed  and  have 
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some  interlocking  of  the  grains  that  account  for  the  higher 
strength . 

The  tangential  stress  at  failure  is  Known  and,  with  the 
friction  angle  determined,  the  radial  stress  at  failure  may 
be  calculated.  It  is  thus  possible  to  plot  these  results  on 
a  p-q  diagram.  The  results  of  Figure  5.4  are  reproduced  in 
Figure  6.16  with  the  hollow  cylinder  data  included.  It  can 
be  seen  from  this  plot  that  the  hollow  cylinder  oil  sand 
strengths  lie  much  closer  to  the  results  for  a  lean  oil  sand 
than  do  the  triaxial  test  data  for  Series  B  and  C. 

6.4.2  Wall  Movements 

The  predicted  wall  movements  result  in  a  good 
comparison  with  the  measured  movements.  The  predicted 
internal  movements  are  lower  than  the  actual  movements  soon 
after  the  commencement  of  yielding.  This  is  due  to  the 
assumption  made  regarding  the  average  radial  stress  in  the 
elastic  zone.  Once  yielding  occurs  and  the  plastic  zone 
begins  to  develop,  the  average  radial  stress  in  the  elastic 
zone  increases  rapidly  and  thus  results  in  an  increase  in 
Young's  modulus.  This  subsequently  results  in  less  elastic 
straining  in  the  plastic  zone. 

The  predicted  movements  are  larger  than  the  actual 
movements  when  the  cylinder  approaches  the  collapse  state, 
this  is  due  to  the  problem  of  modelling  post-failure  strains 
as  well  as  to  the  influence  of  the  overprediction  of  the 
external  wall  movements. 
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The  external  wall  movements  are  slightly  underpredicted 
by  the  model  in  the  elastic  range  and  overpredi cted  as  the 
sample  nears  the  collapse  state.  This  is  a  result  of  the 
determination  of  E  in  the  elastic  zone  and  to  the  rapidly 
increasing  stresses  at  the  external  wall  as  the  plastic 
radius  approaches  the  value  of  the  external  radius. 

The  overprediction  of  the  external  wall  movements  will 
result  in  an  error  in  the  determination  of  the  internal  wall 
movements.  The  sample  radii  are  modified  after  each  stress 
increment  and  the  stress  calculations  are  based  on  these 
modified  radii.  Thus,  as  the  sample  nears  collapse,  the 
stress  distribution  will  be  affected  by  the  error  in  the 
external  boundary  dimension.  The  stresses  will  be  higher 
than  actually  exist,  and  the  internal  wall  movements  will  be 
overpredicted . 

The  importance  of  the  non-linearity  of  the 
stress-strain  curve  is  demonstrated  in  the  results  of  Test 
#A-2.  In  Figure  6.9,  the  predicted  closures  are  far  less 
than  the  measured  ones,  even  in  the  elastic  range  of  the 
material.  When  the  elastic  moduli  are  adjusted  to  account 
for  variation  in  E  at  the  low  stress  level,  the  predicted 
closures  are  much  better  (Figure  6.10).  The  plastic  strain 
increments  are  a  function  of  the  total  elastic  strains  in 
the  plastic  zone,  and  hence  it  is  important  to  correctly 
predict  the  elastic  strains. 

The  incongruity  of  the  stress-strain  curve  for  Test 
#A-2  in  the  low  stress  range  is  difficult  to  account  for. 
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The  loading  cap  should  be  well  seated  after  the 
consolidation  test,  thus  it  does  not  seem  feasible  that  this 
incongruity  is  a  result  of  seating.  The  measured  closures 
also  respond  in  a  similar  manner,  as  can  be  seen  from  Figure 
6.9. 

The  difference  in  closure  curves  between  the  oil  sand 
samples  and  Ottawa  sand  sample  should  be  noted.  The  Ottawa 
sand  responds  in  a  more  ductile  manner  than  the  oil  sand. 

The  closure  curve  for  Ottawa  sand  shows  a  more  non-linear 
behaviour  and  substantial  "flow"  occurs  prior  to  collapse. 

In  the  oil  sand  samples,  the  closure  is  more  linear  and 
collapse  occurs  suddenly.  This  difference  in  behaviour  is 
due  to  the  difference  in  the  shape  of  the  grains  as  well  as 
the  interlocking  of  the  oil  sand  grains. 

The  failure  mode  of  the  hollow  cylinder  samples  has 
been  discussed  in  the  previous  chapter.  However,  the 
post -col  1  apse  mechanism  of  Test  #A-2  has  not  been  discussed. 
Collapse  of  a  cylinder  is  defined  as  a  very  large  volume 
change  of  the  bore  with  a  small  decrease  in  internal 
pressure.  In  this  test,  the  volume  changes  were  allowed  to 
continue  without  shutting  down  the  pressure  systems  (as  is 
the  usual  practise).  After  a  short  time,  the  volume  changes 
stopped  and  further  reduction  of  the  internal  pressure  did 
not  result  in  any  substantial  volume  changes  of  the  bore  or 
the  pore  fluid. 

The  internal  pressure  was  continually  decreased  until 
the  effective  radial  stress  at  the  inside  surface  was  zero. 
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Even  under  this  condition,  volume  changes  did  not  occur. 

The  sample  thus  stabilized  (probably  by  arching)  under  zero 
effective  internal  pressure  and  an  external  pressure  of  10 
MPa. 

Physical  examination  of  the  sample  after  the  test 
showed  a  distinct  failure  surface  typical  for  this  type  of 
test.  The  internal  diameter  was  uniform  except  in  the  zone 
of  the  failure  surface. 

This  arching  phenomenon  needs  to  be  examined 
experimentally  in  further  detail  before  any  conclusions  may 
be  drawn.  The  arching  may  be  a  function  of  the  wall 
thickness  and  is  probably  not  to  be  expected  in  the 
deformation  around  a  shaft  or  tunnel  in  cohesionless 
material.  It  should  be  noted  that  the  shape  of  the  bore  of 
the  Ottawa  sand  sample  after  the  test  was  not 
distinguishable  as  the  sample  had  collapsed  completely. 

Hence  this  arching  phenomenon  is  probably  greatly  influenced 
by  the  shape  of  the  grains  and  the  soil  structure. 


6.5  Conclusions 

The  numerical  model  developed  during  this  research  has 
the  capability  of  predicting  the  deformations  that  occur  in 
a  hollow  cylinder  test  with  reasonable  accuracy.  The 
dilation  parameter  selected  for  the  analysis  was  unity,  the 
case  of  no  plastic  volume  change.  Any  other  values  of  this 
parameter  will  result  in  larger  overpredict  ions  of  the 
deformations  as  the  cylinder  nears  the  collapse  state. 
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The  use  of  a  simplified  analysis  to  determine  the 
stresses  in  the  sample  is  valid  only  when  studying  the 
elastic  deformation  of  the  hollow  cylinders.  Two  different 
assumptions  were  used  in  the  analyses  and  both  gave  similar 
results  in  the  elastic  range. 

These  simplified  analyses  are  not  valid  when 
determining  the  shear  strength  of  the  hollow  cylinders.  The 
shear  strength  as  determined  from  Wu  et  a/. (1963) 
corresponds  with  that  determined  from  the  numerical  model. 
The  model  can  thus  be  used  as  an  analytical  tool  in  the 
reduction  of  the  hollow  cylinder  test  data. 

The  arching  phenomenon  observed  in  Test  #A-2  needs  to 
be  examined  in  further  detail.  This  phenomenon  was  not 
evident  in  the  Ottawa  sand  sample  and  is  probably  a  function 
of  the  grain  shape  and  soil  structure.  The  influence  of 
wall  thickness  has  a  large  effect  on  the  arching  mechanism. 


CHAPTER  7 


CONCLUSIONS 

The  research  herein  has  examined  the  deformation  behaviour 
of  oil  sands.  An  experimental  program  was  conducted  on 
intact  specimens  of  oil  sand.  The  problems  associated  with 
sample  preparation  and  disturbance  have  been  discussed  and  a 
detailed  analysis  of  the  i nterpretat ion  of  hollow  cylinder 
test  data  has  been  made. 

The  research  has  resulted  in  a  better  understanding  of 
the  problems  associated  with  the  testing  of  oil  sands,  as 
well  as  the  deformat iona 1  behaviour  of  this  material.  The 
following  sections  summarize  the  major  findings  of  this 
research . 


7.1  Sample  Preparation  and  Laboratory  Equipment 

The  influence  of  sample  preparation  techniques  on  the 
amount  of  sample  disturbance  that  occurs  has  been  shown.  It 
is  necessary  to  improve  these  techniques  if  reliable  data  is 
to  be  obtained  from  laboratory  testing. 

The  disturbance  that  can  occur  in  the  preparation  of 
samples  is  due  to  the  heat  generated  by  the  coring  and 
trimming  procedures.  The  effects  of  disturbance  can  be 
reduced  by  using  larger  size  samples  to  dissipate  the  heat. 
However,  it  is  necessary  to  improve  these  techniques  so  as 
to  minimize  the  heat  generation. 

Many  of  the  problems  encountered  in  conducting  high 


135 


136 


pressure  tests  have  been  described.  The  hollow  cylinder 
triaxial  device  has  been  used  with  limited  success  and  many 
more  successful  tests  are  required  before  the  equipment  may 
be  called  operational.  The  problem  of  leaks  developing  in 
the  internal  membrane  may  be  minimized  by  using  flexible 
Neoprene  rubber  rings  at  all  interfaces  between  the  sample 
and  the  equipment  (loading  cap,  porous  stones). 

The  high  pressure  triaxial  cell  for  37.5  mm  diameter 
samples  operated  without  any  problems.  This  will  be  an 
extremely  valuable  piece  of  equipment  once  the  sample 
disturbance  problem  for  small  samples  is  solved. 


7.2  Material  Behaviour 

The  laboratory  results  have  shown  that  the  deformation 
behaviour  of  oil  sand  is  path  dependent.  This  path 
dependency  needs  to  be  examined  more  thoroughly.  The 
influence  of  sample  disturbance  will  be  reflected  in  the 
test  results  and  may  not  be  indicative  of  the  in  situ 
behaviour  of  the  material. 

The  oil  sands  deform  in  a  non-linear  manner.  This 
non-linearity  is  the  result  of  irreversible  strains  that 
occur  in  particulate  media  and  is  a  function  of  the  amount 
of  slip  that  occurs  between  particles.  In  the  case  of  the 
hollow  cylinder  tests,  the  reduction  in  the  confinement  of 
the  sample  will  also  result  in  a  non-linear  response. 

The  oil  sands  dilate  when  sheared.  This  dilatancy  is  a 
function  of  the  confining  pressure  and  can  be  suppressed  at 


:  *! 


137 


high  pressures.  This  change  in  mode  of  failure  has  been 
shown  to  be  path  dependent. 

Crushing  of  the  sand  grains  occurs  at  these  high 
pressures  as  indicated  by  the  increase  in  the  amount  of 
fines  after  testing.  Some  crushing  also  occurs  during  the 
consolidation  of  the  samples,  but  these  two  phenomena  have 
not  been  examined  separately. 

An  arching  phenomenon  has  been  observed  in  a  hollow 
cylinder  sample.  Insufficient  test  data  have  made  it 
impossible  to  analyse  this  phenomenon  in  detail,  but  it  is 
dependent  on  the  shape  of  the  sand  grains  as  well  as  the 
interlocking  nature  of  these  materials. 


7.3  Hollow  Cylinder  Test  Data 

The  hollow  cylinder  test  is  ideal  for  examining  the 
deformations  that  occur  around  shafts  and  tunnels.  The 
stress  path  followed  in  this  test  is  similar  to  actual  field 
behaviour,  and  the  influence  of  reduced  confinement  on  the 
deformation  moduli  can  be  readily  examined. 

A  numerical  model  has  been  used  to  show  how  this 
influence  of  confinement  may  be  i ncorporated .  This  approach 
is  simpler  than  the  more  conventional  ones  used  in 
non-linear  analyses.  This  approach  makes  the  hollow 
cylinder  test  a  powerful  tool  for  obtaining  design 
parameters . 

The  i nterpretat ion  of  data  from  a  hollow  cylinder  test 
has  been  examined.  The  assumption  of  axi -symmetr i c 
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deformation  is  valid  as  is  evident  from  the  examination  of 
the  cylinders  after  the  test.  End  restraints  have  no 
significant  influence  on  the  deformed  shape  of  the 
cylinders,  though  the  shape  of  the  failure  surface  is 
influenced  by  these  restraints.  The  uniform  deformation  can 
thus  be  reliably  measured  on  a  volume  change  indicator. 

The  radial  stress  distribution  within  a  hollow  cylinder 
may  be  approximated  as  varying  linearly.  The  stress-strain 
curve  will  give  reliable  deformation  moduli  within  the 
elastic  portion  of  the  curve.  The  shear  strength  of  the 
sample  must  be  determined  from  an  analytical  solution  such 
as  that  given  by  Wu  et  a/.  (1963)  or  from  the  model 
developed  in  Chapter  4.  The  assumptions  made  regarding  the 
radial  stress  distribution  will  result  in  an  incorrect 
determination  of  the  shear  strength. 


7.4  Concluding  Remarks 

The  results  of  this  laboratory  program  have  revealed 
several  traits  of  oil  sand  behaviour  which  should  assist  in 
a  better  understanding  of  this  material.  The  detrimental 
effect  of  sample  disturbance  must  be  overcome  before 
reliable  data  can  be  obtained. 

Many  other  areas  of  research  must  still  be  addressed. 
The  most  important  of  these,  as  far  as  in  situ  recovery 
procedures  are  concerned,  are  an  examination  of  the 
undrained  behaviour  of  oil  sand  and  the  influence  of  high 
temperatures  on  the  mechanical  properties. 
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E last ic  solut ion  for  the  plane-strain  problem 
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The  continuum  mechanics  relationships  for  the  plane-strain 
problem  can  be  found  in  several  elasticity  text  books,  and 
are  summarized  below  in  terms  of  polar  co-ordinates  : 

Equilibrium  Equations: 

+  Fr  =0 

....  A 

=  0 

Hookes  Law: 


€.  =  ^[0,-(7^7)<Tr] 

Strain-Displacement  Equations: 


-V  1^1  +  cht*  _  “ a. 

Yr»  ~  r  00  dr  r 
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Stress  components  in  terms  of  the  Airy  stress  function : 


dr 


=  i0$  ,  i  6^ 


cre  - 


r  0r 


1  0e* 


A  4 

•  •  •  ■  •  r 


r  -  i  dQ  i  0,(fr  _  0/i  0(K 

r#'  ri09  r  0r0e  "  /=ir  V  r  a.  ) 
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The  bi -harmonic  equation  is  derived  by  considering  the 
equilibrium  equations,  Airy  stress  function,  compatibility 
requirements  and  Hookes  Law  : 


V4< t»  =0  ...  .A. 5 

For  the  axi -symmetr ic  solution, -A-  =0.  The  general 
solution  of  the  bi -harmonic  equation  for  this  case  gives  the 
Airy  stress  function  and  hence  the  stress  components  : 


CTr  =  A(1  +2lnr)  +  2  B  + 
C7«=A(2lnr+3)i-2B  — 


.  .  .  .  A  . 


6 


In  the  hollow  cylinder,  the  region  occupied  by  the 
cylinder  is  not  simply  connected  and  hence  the  compatibility 
equations  are  not  sufficient  to  guarantee  single  valued 
displacements.  It  is  therefore  necessary  to  consider  the 
strain-displacement  equations  as  well.  Integration  of  these 
equations  give  radial  displacements  which  must  be 
equivalent.  For  this  condition  to  be  satisfied,  the 


149 


constant  A  in  Eqn.A.6  must  be  zero.  Thus  the  stress 
distribution  in  the  hollow  cylinder  is  given  by  : 

af  -  2B  + 

cre  =  2b  - 

Tr,  *  0 


The  constants  B  and  C  are  solved  from  the  boundary 
condi t ions . 

Boundary  Condi t ions 

The  sign  convention  used  in  the  derivation  is  that 
compressive  strains  and  stresses  are  positive. 

The  boundary  conditions  for  the  hollow  cylinder  test  are  : 


At  r  =  a  ,  G'r  =  P| 


At  r  =  b  ,  CJr  -  Pe 
Pj  = 2B  + 

Pe-  2B+ 


.A. 8 


Elimination  of  2B  from  these  equations  solves  for  C 


C  = 


-aVf^-p^ 


.  .A. 9 


.«  * 
b  -a 
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Elimination  of  C  from  eqn.A.8  solves  for  2B  : 


-pa+p  b* 

2Bn-J - 

ba-  a*  _ A  .  1 0 

Substitution  of  eqns . A . 9  and  A. 10  into  eqn.A.7  gives  the 
solution  for  the  stress  distribution  in  the  hollow  cylinder 


a*  b*(P,  -  P; ) 

r*  ( b*-  a* ) 


oVtP^-P;) 

r‘  (  b*-  a* ) 


=  0 


.  A.  1  1 
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E lasto-plast ic  solut ion  for  the  stress  distr ibut ion  in  the 


hollow  cv 1 i nder  test 
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The  derivations  which  follow  are  based  on  the  following 
assumptions  : 

1.  Linear  Mohr-Coulomb  failure  criterion 

2.  E last ic-perfect ly  plastic  behaviour 

Stress  Distribution  in  the  Plastic  Zone 

The  Mohr-Coulomb  failure  criterion  relates  the  maximum 
and  minimum  principle  stresses  at  failure  as  follows  : 


( 7, »  -  N(7r  +  2cN1/2 


.  .  .  .  B  .  1 


The  equilibrium  equation  for  the  axi -symmetr i c  case  is  : 

-  T-ICJ.-dr) 

dr  r 

•  •  • 

Substituting  Eqn.B.1  into  Eqn.B.2  yields  : 


Integration  of  the  above  equation  determines  the  radial 
stress  distribution  in  the  hollow  cylinder  : 


.  .  .  .  B  .  4 


Boundary  Conditions 

The  boundary  conditions  after  the  initiation  of  failure 
are  given  by  : 


At  f  s  Q  *  “  Pj 


At  r:b,  Cr  *  Pe 


B .  5 


At  r  :  R  ,  (Jr  *  Op 


•  tii 


• ' 


Substitution  of  the  applicable  boundary  condition  into 
Eqn.B.4  solves  for  the  constant  A  : 
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A  =  -! 


P:  (N-1 )  -*-2c  N 


1/2 


.  .  .  .  B  .  6 


N-1 


The  radial  stress  distribution  in  the  plastic  zone  is  thus  : 


<J?  =  j^tP,(N-1)  +  2c  N  ] 


1/2  l r  - 1_(2cN1/2]  •  •  •  •  B  .  7 


N-1 


The  tangential  stress  distribution  is  calculated  from  from 
the  failure  criterion  (Eqn.B.1)  : 


£7,P  =  NCJrP  +  2cN1/2 


[  P.  N  +  2cn3/21 

f  -f"1 

Li  N-1  J 

L  a  J 

2  c  N 

N-1 


,1/2 


.  .  .  .  B  .  8 


The  radial  stress  at  the  elastic-plastic  boundary  is 
calculated  from  Eqn.B.7  and  the  boundary  condition  (Crrs(Tp 
at  r  =  R  )  : 


Stress  Distribution  in  the  Elastic  Zone 
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The  solution  is  similar  to  that  derived  in  Appendix  A 
(Eqn.A.11).  However,  the  boundary  conditions  are  modified 
as  fol lows  : 


At  r  =  R  ,  CTr  =  Op 


At  r  =b  ,  CJr  -  Pe 

Thus  in  Eqn.A.11,  if  'a'  is  replaced  by  '  R  '  and  '  p.  '  by 
'  ,  the  solution  is  given  by  : 


-CTpR* 

b‘-  R2 


Rzb2(pp-(JD) 
rx(  b*-R£) 


R2b2(Pe-ap) 

+  r1  (  b2-  Rx) 


.  B  .  1  0 


=  0 


The  stresses  in  the  elastic  zone  cannot  be  determined  until 
the  radius  of  the  plastic  zone  is  found. 


Radius  of  the  Plastic  Zone 

It  is  assumed  that  the  radial  stresses  are  continuous 
at  the  elastic-plastic  boundary.  From  the  elastic  solution 
it  can  be  shown  that  : 


b*  -  Pjg2 
b1-  <r 


)  =  2C70 


•  •  • 


.B.  1  1 
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This  equation  is  valid  for  all  values  of  the  radius  in  the 
elastic  zone.  Thus  : 


(  (Jr  ) r  -  p  =  Op  -  2  Oq  Of 


.  B  .  1 2 


At  the  plastic  boundary,  failure  has  just  been  reached  and 
hence  the  failure  criterion  can  be  applied  : 


*  V  s  R  *  0QR  =  N  0*p  +  2  c  N1/2 


rr  J/2 

UaR  -  2  c  N 

N 


.  B  .  1  3 


Elimination  of  from  eqns.B.12  and  B.13  : 


2  CT0  -  2c  N 

N  +1 


1/2 


. . . . B .  14 


The  radial  stress  at  the  boundary  from  the  plastic  solution 
i  s  : 


.B.  15 


Equating  Eqns.B.14  and  B.15  solves  for  the  radius  of  the 
plastic  zone  : 


R 


a 


N*1 


[ 


1/2 

cr0(N-n  +  2c  n 

Pj(N-l)  +  2cN1/2 


V 


.  B  .  1 6 
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Average  Radial  Stress  in  the  Cylinder 

The  total  radial  force  in  the  cylinder 


is  given  by  : 


/ 


.  B  .  1 7 


Now , 


for  a  cohesionless  material 


r  RN-  aN 

L  NaN-’ 


] 


b*  -  Op  R1 


bV(P,-(JD) 

.  r*  (  bx-  Rl) 


R 


Peb(b-R)  -  Op  R  (R-b) 
b+  R 

Substituting  for  (Jp  (Eqn.B.9)  and  rearranging  : 


Fr 


a.  (b-R)  RN  1 
N  (b+R)  '  „N-1  J 


Peb 


(fa-R)  _ B .  18 

(b+R) 


The  average  radial  stress  in  the  cylinder  is  then  defined  as 


Jll. 

b  -  a 


^r  lave 


•  iii 


B.  19 


1  li 
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A  Versat i 1e  Hollow  Cv 1 i nder  Tr iaxia 1  Dev i ce  (a  copy  of  a 

paper  by  M.B.  Dusseault,  1981) 
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by 
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ABSTRACT 


A  high-pressure,  thick-walled  hollow  cylinder  test  device  has 
been  perfected  to  explore  the  consequences  of  stress  paths  which  are 
typical  of  pressuremeter  configurations,  tunnelling,  or  shafting  in 
soils  materials.  Several  features  have  been  incorporated  to  facilitate 
data  gathering,  and  the  device  is  readily  adaptable  to  Kq  or  standard 
compression  and  extension  triaxial  tests. 

The  device  has  been  used  on  tests  in  artificial  dense  sand 
specimens  constructed  by  a  freezing  technique,  and  will  be  used  in  an 
extensive  study  of  the  behavior  of  oil  sands.  Typical  results  of  tests 
on  dense  Ottawa  Sand  (C-109)  are  presented. 
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INTRODUCTION 

Greater  design  demands  and  the  increasing  capabilities  of 
numerical- model!  ing  techniques  in  geomechanics  have  limited  the  useful¬ 
ness  of  the  more  simplistic  constitutive  laws  of  soil  behavior.  There 
is  a  need  for  laboratory  testing  techniques  which  can  follow  stress 
paths  similar  to  those  which  soils  undergo  in  nature.  The  popularity  of 
tests  such  as  the  direct  shear  and  triaxial  compression  tests  lies  in 
their  simplicity  and  the  widespread  need  for  strength  data  for  limit 
equilibrium  analysis;  however,  these  tests  do  not  provide  data  adequate 
for  all  deformation  behavior  analyses. 

The  hollow  cylinder  triaxial  device  (HCTD)  has  been  designed 
and  built  to  provide  data  from  tests  with  careful  stress  path  control. 
The  HCTD  can  also  serve  as  a  model  of  shafting  in  a  uniform  horizontal 
stress  field,  to  explore  material  behavior  under  various  radial  stress 
gradients,  and  (with  minor  modifications)  for  a  variety  of  standard 
tests. 

THE  UNIVERSITY  OF  ALBERTA  HIGH  PRESSURE  HOLLOW  CYLINDER  TRIAXIAL  DEVICE 

The  device  was  designed  in  response  to  the  needs  of  oil  sands 
research  and  development.  Mine-assisted  in  situ  projects  will  require 
shafts  and  tunnels  to  gain  access  to  deep  deposits  (Dusseault,  1978), 
and  massive  surface  injection/production  jm  situ  projects  require 
knowledge  of  complex  material  behavior  (Imperial  Oil  Limited,  1978)  such 
as  fracture  toughness. 

The  initial  step  in  the  development  of  the  HCTD  was  to  build  a 
special  top  cap  for  a  standard  triaxial  device  (Figure  1).  The  top-hat 
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internal  membrane  is  formed  by  molding  and  curing  latex  on  an  aluminum 
mold,  and  the  Teflon  membrane  seat  prevents  internal  membrane  rupture. 
This  low-pressure  device  (1.7  MPa)  does  not  permit  axial  stress  control 
or  true  plane  strain,  therefore  a  true  HCTD  was  constructed.  The  first 
HCTD  (which  is  very  similar  to  the  second  described  more  fully  below) 
has  a  maximum  confining  pressure  (pQ)  of  6.3  MPa,  and  has  a  top  cap 
design  (Figure  2)  which  permits  axial  fluid  flow  but  provides  poor 
sealing  capability  at  pressures  greater  than  5.5  MPa.  This  first  device 
has  proven  very  useful  in  the  evaluation  of  the  behavior  of  dense  sands 
and  is  still  in  use;  however,  because  the  oil  sand  displays  high 
strengths  and  a  curvilinear  Mohr-Coulomb  failure  envelope  (Dusseault  and 
Morgenstern,  1979),  higher  outside  pressures  are  required  to  carry  a 
hollow-cylinder  test  to  the  state  of  plane  strain  cylinder  collapse. 

The  second  HCTD  (Figures  3  and  4)  is  a  high-pressure  thick- 
walled  cylinder  testing  apparatus  with  the  following  capabilities: 

(1)  high  pressure  range  (pQ  =  20  MPa  with  FS  of  about  1.5); 

(2)  sample  size  of  101.6  mm  outside  diameter,  50.8  run  inside 
diameter,  200  to  240  mm  in  length  (depending  on  membranes 
available) ; 

(3)  independent  pressure  control  and  volume  change  measurements  on 
inside,  outside  and  pore  liquids;  and 

(4)  a  101.6  mm  ram  with  a  low-friction  seal  and  Teflon  ram  guide 
to  allow  plane  axial  strain  with  passive  stress  measurement, 
independent  axial  stress  control  (K  control),  or  specified 
strain  rate  with  stress  measurement. 
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To  increase  the  versatility  of  the  HCTD  the  following  features 
were  included  in  the  design: 

(1)  '  the  apparatus  is  built  of  standard  available  corrosionless 

material  and  0-ring  seals; 

(2)  the  cell  base,  body,  and  top  can  be  used  directly  for  other 
sample  inside  and  outside  diameters  which  require  only  minor 
redesign  (cap  and  pedestal,  ram  guide,  ram,  or  membranes); 

(3)  a  ram  guide  and  ram  for  easy  conversion  to  conventional  tri- 
axial  compression  tests  is  available;  and 

(4)  sufficient  interior  clearance  has  been  provided  for  heating 
and  cooling  coils,  or  for  mounting  of  diametral  strain 
measurement  devices. 

This  HCTD  is  equipped  with  17.5  MPa  working  pressure  no¬ 
volume-change  valves  and  3.175  mm  brass  tubing;  a  44.5  kN  axial  load 
cell  (az);  a  14  MPa  fluid  pressure  transducer;  two  high-pressure  (14 
MPa)  systems  with  two  in-line  14  MPa  volume  change  indicators  allowing 
pressure  control  of  pQ,  p^,  or  monitoring  of  A V  ,  AV^;  and  two  low- 
pressure  regulator-controlled  pressure  systems  with  individual  volume 
change  indicators  for  the  pore  fluid  or  for  lower  pressures  (<  1.4  MPa) 
on  p..  or  pQ.  For  testing,  the  HCTD  is  mounted  in  an  extremely  stiff 
plane  strain  test  frame,  or  in  a  90  kN  capacity  infinitely  variable 
strain  rate  test  frame.  The  cell  has  been  tested  to  17  MPa  and  has  been 
used  in  tests  on  dense  sands  to  confining  stresses  (pQ)  of  12.6  MPa. 
Tests  with  pQ  less  than  about  6  MPa  use  latex  exterior  membranes  (1.5  mm 
thick)  molded  in  our  laboratory;  higher  pressures  require  Neoprene 
exterior  membranes  (3  mm  thick).  High-pressure  membranes  require 
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special  molds  and  high-pressure  and  high-temperature  injection  of  the 
Neoprene.  The  molds  have  been  made  in  the  University  of  Alberta  and 
Neoprene  injection  is  performed  commercially. 

The  fully  assembled  fluid-filled  HCTD  is  heavy  and  requires  a 
small  crane  for  handling.  Individual  components  can  be  lifted  by  one 
person,  but  assembly  is  simplified  if  a  second  person  is  available  for 
about  15  minutes. 

SAMPLE  ASSEMBLY  AND  MOUNTING 

Specimens  of  significant  cohesive  strength  can  be  prepared  and 
mounted  directly  in  the  HCTD,  but  fabricated  specimens  of  dense  sand  or 
oil  sand  specimens  require  special  preparation  techniques. 

Ordinary  specimens  are  trimmed  to  101.6  mm  external  diameter, 
placed  in  a  split  retaining  cylinder  and  cored  with  a  50.8  mm  diamond 
coring  device.  Oil  sand  specimens  must  be  trimmed  to  diameter  in  a 
lathe  by  use  of  tungsten-carbide  tipped  bits  (Dusseault  and  Morgenstern, 
1978)  in  a  cold  room  at  temperatures  of  -25°C  to  prevent  sample  dete¬ 
rioration.  The  core  drilling  must  be  carried  out  with  the  oil  sand  at 
very  low  temperatures  and  heat  from  the  coring  process  must  not  be 
allowed  to  build  up. 

Because  of  the  top  hat  design  of  the  inner  and  outer  mem¬ 
branes,  a  special  jig  (Figure  5)  has  been  built  to  allow  thick-walled 
sand  cylinders  to  be  assembled  under  water  directly  on  a  vibratory  (60 
Hz)  table.  Shock  loading,  static  top  loads,  and  rodding  during  assembly 
are  aids  to  densif ication.  A  porous  stone  is  placed  on  the  top  of  the 
specimen,  and  a  perforated  lucite  top  disc  is  placed  on  the  stone.  The 
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bottom  of  the  jig  is  placed  in  a  cold  bath  (-18°C)  in  a  cold  room 
(-2°C)  and  a  static  load  of  250  N  is  placed  on  the  Lucite  top  disc'. 
One-dimensional  freezing  takes  place  with  excess  fluid  expelled  through 
the  top.  When  frozen,  the  jig  is  disassembled  (special  techniques  are 
required)  and  the  sample  is  mounted  in  the  HCTD  in  the  cold  room. 

Special  methods  (e.g.  back-pressures  and  pore  suction)  are  necessary 
during  initial  flooding  of  the  cell  and  the  internal  membrane  to  prevent 
collapse  of  the  dense  sand  cylinder. 

Many  small  steps  have  been  omitted  from  this  description.  A 
comprehensive  paper  on  sample  and  test  techniques  is  being  prepared. 

CAPABILITIES  OF  THE  UNIVERSITY  OF  ALBERTA  HCTD 

A  variety  of  test  types  can  be  pursued  using  the  HCTD.  The 
major  stress  paths  and  the  ones  of  greatest  interest  deserve  particular 
mention. 

Radial  Strain  Tests 

For  plane  strain  cylinder  collapse  (PSCC),  three  basic  con¬ 
figurations  are  possible;  p^ ,  pQ  or  u  may  be  monotonical ly  changed  until 
failure  occurs. 

0)  pQ  increasing;  p^,  u  constant.  A  sample  is  stabilized  at  the 

stress  conditions  pQ  =  p^  >  u,  then  pQ  is  increased  to  failure 

AV-,  A V  ,  e  ,  a  are  measured. 

1  u  z  z 

(2)  p.  decreasing;  pQ,  u  constant.  A  sample  is  stabilized  at  the 
stress  conditions  pQ  =  p^  »  u,  then  p^  is  decreased  to 
failure;  AVq,  AVu>  ez,  cz  are  measured. 
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(3)  u  increasing;  pQ,  p..  constant.  A  sample  is  stabilized  at  the 

stress  conditions  P0  >  P-j  >  u>  then  pQ  and  p..  are  held 

constant  while  u  is  increased  to  failure;  AV  ,  A V.,  e  ,  a  are 

o  1  z  z 

measured. 

These  three  test  configurations  are  statically  determinate 
(Timoshenko  and  Goodier,  1951),  and  therefore  stresses  may  be  calculated 
directly  by  specifying  a  failure  criterion  only  (Wu  et  aK,  1963). 

For  plane  strain  cylinder  burst  (PSCB)  conditions,  statically 
determinate  tests  may  be  performed  in  the  HCTD  which  are  directly 
analagous  to  the  PSCC  tests.  Note  that  u  can  at  no  time  be  greater  than 
pQ  or  p^,  but  that  p^  is  greater  than  pQ  during  PSCB  conditions. 

Because  failure  is  catastrophic  in  PSCC  and  PSCB  tests,  constant  vigi¬ 
lance  is  required  near  failure  to  shut  down  pressure  lines.  All  tests 
can  be  performed  in  an  undrained  mode  with  measurement  of  pore  pressure 
rather  than  pore  volume  change. 

Axial  Strain-Controlled  and  Stress-Controlled  Tests 

Test  configurations  where  failure  is  approached  by  straining 
in  an  axial  direction  are  in  general  not  statically  determinate  (Wu  et 
al . ,  1963).  Since  pQ,  p^,  and  u  may  be  kept  constant,  complete  volume 
change  data  may  be  collected  during  testing.  Normal  stress  gradients 
(p  >  p^)  or  reverse  gradients  (p^  >  pQ)  may  be  used  across  the  radius 
of  the  specimen.  Because  the  ram  is  the  same  diameter  as  the  specimen, 
compression  or  extension  tests  may  be  performed. 

The  introduction  of  a  servo-mechanism  or  a  constant-stress 
load  device  along  the  specimen  axis  will  permit  stress  control  or  creep 
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testing  of  thick-walled  hollow  cylinders.  A  stress-control  feature  is 
necessary  if  the  device  is  to  be  used  as  a  shaft-modelling  apparatus*. 


Compressibility  Tests 

Because  the  ram  is  the  same  diameter  as  the  cylindrical  speci¬ 
men,  a  no-axial-strain  compressibility  test  can  be  performed.  The  basic 
triaxial  cell  yields  an  all-round  compressibil ity  (Aa-j  =  Aa^  =  Aa3);  the 
oedometer  provides  the  condition  A az  =  Aa^  >  (Aa^  =  A a^)  during  loading; 
and  the  HCTD  provides  the  condition  (Aa^  =  Aa^)  >  (A =  A a  )  during 
loading  with  no  axial  yield  providing  that  P-j  =  PQ  >  u  for  all  stress 
increments.  For  the  plane  axial  strain  HCTD  compressibility  test,  AVu, 
ez  and  az  are  measured  after  each  stress  step.  Compressibility  tests 
may  also  be  performed  at  a  constant  radial  stress  ratio;  i.e.  (p  -  u )/ 
(p..  -  u)  as  a  positive  constant  less  than  the  failure  ratio.  More 
complex  stress  paths  are  readily  envisioned,' 


TYPICAL  TEST  DATA 

Cyclic  plane  axial  strain  compressibility  tests  v/ere  performed 
on  a  dense  Ottawa  Sand.  To  confirm  the  assumption  of  plane  strain, 
axial  displacement  was  measured  as  well  as  axial  load.  A  constant  pore 
pressure  permitted  accurate  determination  of  pore  fluid  volume  changes, 
and  pQ  was  identical  to  p^  for  all  stress  increments.  The  data  (Figure 
6)  indicate  small  axial  strains  (some  compliance  must  exist  because 
of  the  load  cell  and  confinement  frame);  therefore,  plane  strain  con- 
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ditions  are  approached.  The  coefficients  of  volume  compressibility  are 
presented  in  Table  1 . 

Two  types  of  cylinder  collapse  tests  on  dense  Ottawa  Sand  have 
been  performed:  PSCC-Type  1  and  PSCC-Type  2.  Figures  7  and  8  show  data 
generated  during  typical  tests.  Volume  changes  of  pore  fluid  are 
expressed  as  percentages  of  pore  volume  and  may  be  converted  directly  to 
total  volume  change  by  multiplying  by  the  porosity.  These  two  tests 
correspond  to  increasing  J1  or  decreasing  Jl,  and  the  volume  change 
behavior  of  the  specimens  is  therefore  quite  different.  To  ensure  that 
the  specimens  were  very  similar,  cyclic  consolidation  (compressibility 
tests)  to  12.7  MPa  was  performed  on  both  specimens  until  elastic  (recover¬ 
able)  response  was  approximated.  Further  comments  and  interpretation  of 
the  test  data  are  reserved  for  a  future  publication  on  the  strength  and 
volume  change  behavior  of  oil  sand  and  very  dense  Ottawa  Sand. 

CONCLUDING  REMARKS 

The  hollow  cylinder  triaxial  concept  is  not  original  (Kirkpatrick, 
1957;  Haythornthwaite,  1960;  Wu  et  a]_. ,  1963;  Handin  et  al_. ,  1967; 

Jaeger  and  Cook,  1969);  however,  the  University  of  Alberta  high-pressure 
HCTD  incorporates  a  number  of  new  features  permitting  a  great  variety  of 
test  configurations  including  a  simple  conversion  to  standard  triaxial 
testing.  The  HCTD  was  designed  primarily  for  plane  strain  (ez  =  0), 
high-pressure,  thick-walled  cylinder  tests  of  the  collapse  or  burst 
type,  and  for  cyclic  compressibility  tests.  Other  test  configuations 
with  different  materials  are  being  planned. 
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LIST  OF  SYMBOLS  AND  ABBREVIATIONS 


P  »  P- 
Ko 


u 


change  (in  volume,  pressure  or  stress) 
outside  and  inside  pressures 
pore  pressure 


a  ,  a, ,  a?,  mean  axial  stress,  major,  intermediate,  and  minor 
z  principal  stresses 


J1 

K_ 


V  V  vi 


z 

HCTD 

PSCC;  PSCB 


°1  +  °2  +  a3 


-  first  stress  invariant 


principal  stress  ratio 

pore,  outside,  and  inside  volumes 

axial  strain 

hollow  cylinder  triaxial  device 
plane-strain  cylinder  collapse;  cylinder  burst 
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DRAINED  HOLLOW  CYLINDER  COMPRESSIBILITY 
TEST  DATA  (Api=ApQ) 


p0(pi) 

Aa 

z 

Ae 

z 

AV 

u 

mv (=ApAV^ 

kPa 

kPa 

xlO"3 

V 

V  ' 

xlO”3 

kPa  lx  10 

696 

-> 

1390 

+  268 

-0.69 

-3.38 

4.87 

1390 

-> 

2774 

+  584 

-0.37 

-3.22 

2.33 

2774 

-v 

4143 

+  596 

+  0.61 

-2.09 

1.53 

4143 

-> 

5650 

+  673 

+  0.95 

-1.85 

1.23 

5650 

-> 

6926 

+  513 

+1.01 

-1.61 

1.26 

6926 

-> 

8316 

+563 

+1.01 

-1.20 

0.86 

8316 

-* 

6929 

-494 

-0.89 

+  0.72 

0.52 

6929 

4150 

-663 

-2.31 

+  2.09 

0.75 

4150 

-> 

1396 

-1240 

-0.62 

+  4.82 

1.75 

1396 

-v 

4139 

+  939 

+0.41 

-3.86 

1.41 

4139 

-> 

6926 

+  881 

+  2.26 

-2.81 

1.01 

6926 

-> 

8319 

+  543 

+1.20 

-1.21 

0.87 

8319 

-> 

6920 

-431 

-0.89 

+  0.64 

0.52 

6920 

4144 

-620 

-2.34 

+  2.09 

0.75 

4144 

-> 

1392 

-1228 

-0.69 

+  4.90 

1.78 

1392 

-> 

4149 

+  945 

+0.53 

-3.94 

1.43 

4149 

6922 

+  865 

+  2.17 

-2.81 

1.01 

6922 

-V 

8328 

+477 

+  1.31 

-1.13 

0.81 

Note : 


AV 

V 


Since  all  strains  are  known,  and  assuming  elastic 
isotropic  behavior 

Gj+Cj+Cj  and  e  2  =  s3 


2Ap2-ApAa  -A a  2 

E  =  Young's  Modulus  =  — - 

Ap^  -  Aez  (Aaz+2Ap) 


and  v 


Poisson's  Ratio 


Aa 


AV 

z  V 


2Ap 


AV 


-  As  (Aa  +2Ap) 

z  z 

-  2Ae  (A a  +2Ap) 

z  z 


V 
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Figure  1.  Top-Cap  and  Top-Hat  Membrane  for  a 
Conventional  101.6  mm  Triaxial  Cell. 
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Figure  2.  Top-Cap  and  Ram  Assembly  for  6.3  MPa 
Hollow  Cylinder  Triaxial  Apparatus 
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Figure  3.  Exploded  Sectional  View  of  the  High  Pressure  Triaxial  Cell. 
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Figure  4: 


Hollow  Cylinder  Triaxial  Cell  Fully  Instrumented 
in  a  Plane  Strain  Frame 
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Figure  5.  Mold  for  Hollow  Cylinder  Sample  Assembly 
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EFF.  STRESS  KPA 


FIGURE6:  CYCLIC  PLANE  AXIAL  STRAIN  COMPRESSIBILITY  DATA  (PO=PI) 
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FIGURE7:  PLANE  STRAIN  CYLINDER  COLLAPSE  BY  INCREASE  OF  EXTERNAL 
PRESSURE 


APPENDIX  D 


Computat ion  of  the  strains  in  the  hollow  cvl inder  test  (a 


copy  of  pp. 42-46  of  a  thesis  by  D.C.  Procter,  1957) 
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a  2  ”  T  should  vary;  the  size  of  the  variation  will  be  an  indic¬ 
ation  of  the  size  of  th'e  distribution  error# 

b)  ihe  size  of  tne  distribution  error  can  also  be  assessed  by 

comparing  the  peak  strength  values  of  tests  conducted  with  CT~  =  (X0 

with  equivalent  results  from  tests  with  0*  2  =  CT  ,  for  the  same 

^  T 

value  of  N  and  the  same  cize  of  specimen. 

c)  The  overall  error  in  stress  analysis  for  both  sizes  of  cylinder 
can  be  estimated  by  comparing  plane  strain  peak  strengths  from  both 
^2  =  ^  T  ®  2  ~  ®  R  tests  with  the  equivalent  strength  obtained 
in  the  conventional  plane  strain  apparatus. 

Using  an  approach  of  this  nature  it  is  possible  to  obtain  some 
indication  of  the  size  of  the  errors  which  are  likely  to  affect  the 
results. 

Finally  peak  strength  results  obtained  from  thismean  stress  analysis 
are  compared  to  the  same  results  analysed  by  the  method  suggested 
by  Haythomthwaite  (1966).  This  serves  as  a  further  check  on  the 
accuracy  of  the  analysis,  and  is  discussed  in  section  3:4:2. 

30  Mean  Principal  Strain  Commutation 

It  is  necessary  to  compute  from  measured  displacements  and  volume 
changes  the  principal  strain  increments  during  an  imposed  increment 
in  the  azcial  direction.  Since  the  calculated  displacements  can  only 
be  mean  values,  it  is  clear  that  only  mean  strains  can  be  computed 

with  any  accuracy.  It  is  assumed  that  the  principal  stress  and  strain 
directions  are  mutually  coincident  with  the  axes  of  specimen  symm¬ 
etry.  It  is  noted  that  this  computation,  while  being  relatively 
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simple,  is  nevertheless  tedious,  and  consequently  the  following  analysis 
is  presented  in  a  fora  which  can  be  easily  programed  for  a  digital 
computer. 

It  is  required  to  investigate  the  mean  axial,  radial  and  tangential 
principal  strains  8  6  and  £  ^  and  also  the  volumetric  strain  v  , 

by  considering  the  changes  measured  in  the  specimen  length  L,  specimen 
volume  Vs,  and  the  inner  bore  volume  Vb. 

For  the  sake  of  continuity  it  is  helpful  to  state  a  notation  con¬ 
vention: 

Suffix  (  )i  -  indicates  the  value  of  a  quantity  after  the  iL^ 
increment,  i.e.  (L)i. 

*tii 

Suffix  (  )j  indicates  the  value  of  a  quantity  after  the  j 
increment,  i.e,  (L)j,  and  j  =  i  +  1. 

Suffix  (  ) o  —  indicates  the  initial  values  of  a  quantity, 

i.e.  (L)o. 


Prefix  d  -  indicates  the  increment  of  principal  strain  occurring 
during  the  increment,  i.e.  d£^. 

Prefix  D  -  indicates  the  total  change  which  has  occurred  in  a 
quantity  after  a  specified  increment,  i.e.  (DL)^. 

Compressive,  natural  strains  are  considered  positive,  and  the  snail 

strain  approximation  is  assumed  to  apply: 

dv  =  aeA  +  aea  +  deT  I  (300) 

The  mean  radial  strain  must  now  be  defined.  It  is,  however,  more 


convenient  to  define  the  mean  strain  increment  rather  than  the  total 


V 
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-e an  strain,  whence: 


,th  . 


d  €.  -  near,  chnnma  of  specimen  wall  thickness  during  ,j  increment 

mean  wall  thickhess  during  j  increment. 


Hence  if  y  is  the  mean  wall  thickness: 


d6a  = 


(y)j  -  (y).i 

(y)i  +  (y)j 


(3:31) 


Clearly  y  =  b  -  c,  where  b  and  c  are  the  outer  and  inner  radii  of 


the  specimen 


Similarly  d  £  ^  = 


(l)i  -  (L).i 


(3:32) 


U.)i  +  (L>2 

It  now  remains  to  calculate  the  progressive  values  of  L,  b  and  c. 
This  approach  is  a  convenient  one  since  these  values  are  also  required 
for  mean  stress  calculations  (section  3:2). 

If  the  total  volume  of  specimen  and  bore  is  Vz,  where  =  Vb  *  Vs, 
it  follows  that: 


(Lj;  =  (lJo  —  (OL.)  j 


(3:33) 


/ 

l 


fvb]i 
ttU-1  i 


and  similarly  for  (L)j>  (b)j  and  (c)j. 
Hence  from  equations  3:32  and  3:33 


(3:34) 

(3:35) 


Mi  -|Dl)i 
Mo  -(dl)j  -fOL)i 


(3:36) 


' 


. 
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ana  from  equation  30 1: 


de, 


(b)  i  — (b| )  -  f(c)i 

-  &  J 1 

to'  +(»jjj  -(Cq)i 

+  Cc)j]  

X  2 


(3:37) 


where,  from  equation  3: 34 


(b) i  = 

ft)j  = 

and  from  equation  3:35 

(c) i  = 


-  (DVb)  i 


(v s) o  +  (vb)  o  —  (0  Vs)  i 
7T(lL)o  -  (DL)  i) 

(vsjo  +  (vb^s  -  (ovs)j  —  (bvt>)j 

TT((L)o  -  (OQj) 


=  /  f  (y b ;  0  -(DVbji  ~j 

J  L  iT((L)o-(og  i  J 
(c)  j  =  /  r  (vb)o  -Covb)j  -j 

J  L  n(L)0-pL)i)  J 


Also  from  equation  300 


d£  =  dv  -  d  £  ,  -  d  Q  _ 
T  A  R 


where 


dv 


r  (ovs)j  -  (dvs)  <  -i 

i_2(vs)0-(DVsJj-(DVs)i  J 


x  2 


(3:33) 

(3:39) 


Finally  total  strain  calculations  can  be  made  on  a  cumulative  basis: 


.1 


(£  Jj  =  z 
*  1 

(de2) 

(3:40) 

(eA)j  =  i 

A  1 

(eT)j  =  z 

1  1 

(d£A) 

(3:41) 

(de  T) 

(3: 42) 

• 

(v)j  =  z 

•  i 

(dv) 

(3-43) 

From  equations  306  to  3:43  inclusive  all  strain  quantities  required 
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for  analysis  can  be  obtained, 

3:4  Alternative  Failure  Criteria  for  Soil 

3:4*1  Conventional  Approach 

The  philosophy  concerning  the  conception  of  various  failure  criteria 
is  well  known*  and  need  only  be  briefly  dealt  with. 

A  general  stress  state  O’  ^*  O'  O’ ^  may  be  represented  in  3“dimensional 
stress  space  by  two  vector  components*  the  hydrostatic  stress  component 
along  the  hydrostatic  axis  (Q^  =  O' 2  =  ^3)  the  deviatoric  stress 
component*  being  situated  in  the  octahedral  plane  perpendicular  to  the 

hydrostatic  axis  and  containing  the  limit  of  the  hydrostatic  stress 
component  (Fig.  3c) •  The  locus  of  all  such  stress  states  which  are 
situated  in  the  same  octahedral  plane  and  at  which  failure  of  a  partic¬ 
ular  soil  occurs,  is  a  failure  envelope  bounding  all  possible  stress 
states  that  may  occur  for  that  soil  in  that  plane.  Since  the  stress  level 
at  yield  is  dependent  on  the  mean  stress  applied  to  the  soil* . the  size 
of  the  envelope  is  directly  dependent  on  the  hydrostatic  stress  component 
associated  with  the  octahedral  plane  containing  the  envelope.  Consequent¬ 
ly  generating  lines*  meeting  at  the  space  origin  for  the  case  of  cohes¬ 
ionless  materials,  may  be  drawn  through  the  envelopes  on  all  octahedral 
planes  to  form  a  conical  shaped  failure  surface.  These  surfaces  are 
described  in  general  by  a  governing  mathematical  equation  relating 
ffl*  2’  ^3  and  0,  and  it  is  the  form  of  the  equation  which  determines 
the  resulting  shape  of  the  surface. 

The  three  failure  criteria  normally  considered  in  the  case  of  soil 


are  the  Mohr-Coulomb*  the  Extended  Tresca  and  the  Extended  Von  Mises, 


’ 


■ 


APPENDIX  E 


Cal ibrat ion  of  the 


Hoi  low  Cv  1  i rider  Tr  iax ia  1  Dev i ce 
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The  HCTD  has  to  be  calibrated  for  volume  change  correction 
and  ram  friction.  The  volume  changes  that  are  measured  are 
subjected  to  three  sources  of  error  : 

1.  expansion  of  the  cell  and  pressure  lines 

2.  compression  of  the  water  in  the  cell,  lines  and 
measuring  systems 

3.  air  going  into  solution 

The  axial  load  on  the  sample  is  measured  outside  the 
cell,  hence  the  vertical  load  acting  at  the  top  of  the 
sample  has  to  be  corrected  for  the  friction  developed 
between  the  ram  and  the  cell. 

Calibration  of  the  Bore  Volume  Change  Apparatus 

The  sample  used  for  calibration  purposes  is  a  hollow 
cylinder  aluminum  one  of  similar  dimensions  to  that  of  a 
test  specimen.  The  sample  is  assembled  in  an  identical 
manner  as  previously  described  for  the  soil  samples,  with 
the  same  system  of  membranes,  0-rings  and  clamps. 

The  HCTD  was  assembled  in  the  plane-strain  frame  and 
the  bore  pressure  applied  by  means  of  the  nitrogen  pressure 
system.  The  bore  pressure  was  increased  in  increments  of 
approximately  1.4  MPa  and  the  volume  changes  recorded. 

Volume  changes  due  to  the  deformation  of  the  sample  is  small 
and  were  ignored.  The  results  of  the  calibration  are  shown 
in  F i gure  C . 1 . 

Calibration  of  the  Cel  1  Vo  1 ume  Change  Apparatus 


BORE  PRESSURE  (KPfl) 
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TOTAL  VOL.  CHANGE  (CC) 


Figure  E.1  Bore  Volume  Change  Correction  Curve 
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The  cell  pressure  was  applied  in  a  similar  procedure  as 
described  above.  At  the  end  of  each  increment,  the  bore 
pressure  was  increased  to  approximately  the  same  pressure  in 
order  to  minimize  the  sample  deformation.  Volume  changes 
due  to  sample  deformation  were  again  ignored.  The  results 
of  this  calibration  are  shown  in  Figure  C.2. 

Calibration  for  Ram  Friction 

This  calibration  provided  several  difficulties.  Two 
load  cells  would  be  required  to  measure  the  ram  friction, 
one  on  top  of  the  sample  and  the  other  on  the  outside  of  the 
cell.  The  load  cell  inside  the  cell  would  have  had  to  be 
designed  to  withstand  cell  pressures  of  up  to  14  MPa  and 
would  also  have  to  be  specially  designed  as  an  annular  ring 
in  order  to  accomodate  the  internal  membrane  set  up.  It  was 
thus  decided  to  find  an  alternate  method  of  measuring  the 
ram  friction. 

A  machined  aluminum  plate,  4.05  in.  in  diameter,  was 
attached  to  the  loading  cap  and  initially  supported  on  an 
aluminum  block  sitting  on  the  base  of  the  cell.  Upon 
application  of  cell  pressure,  the  loading  cap  and  ram  would 
be  activated  and  an  external  axial  load  could  be  measured  on 
the  load  cell.  This  method  is  equivalent  to  actual  test 
loading  conditions  but  does  not  allow  the  investigation  of 
ram  friction  dependency  on  cell  pressure. 

The  aluminum  plate  is  attached  to  the  loading  cap  by 
clamping  a  small  rubber  membrane  to  the  plate  and  cap  and 


. 


. 


, 


CELL  PRESSURE  (KPfl) 
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TOTRL  VOL.  CHANGE  (CC) 


Figure  E.2  Cell  Volume  Change  Correction  Curve 
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sealing  the  membrane  at  the  ram  by  means  of  0-rings  and 
clamps.  The  cell  pressure  is  then  applied  and  measured  on  a 
transducer.  The  load  cell  is  read  at  each  pressure 
increment.  Results  of  the  test  are  shown  in  Figure  C.3. 

This  calibration  curve  shows  that  ram  friction  is  almost 
independent  of  cell  pressure. 

Correct i on  for  Membr ane  Effects 

The  effect  of  membrane  thickness  on  the  strength  of  the 
samples  has  been  ignored.  Bishop  and  Henkel  (1976)  have 
shown  this  effect  to  be  small  for  1.5  in.  diameter  samples. 
They  also  show  that  for  undrained  tests  no  hoop  tension  is 
induced  in  the  membrane  and  that  the  correction  is  applied 
to  the  axial  stresses  and  not  the  lateral  pressure.  For  4 
in.  diameter  samples,  the  correction  is  strain  dependent  and 
hence  zero  for  plane-strain  conditions  in  the  vertical 
direction  (as  is  the  case  in  this  testing  program) . 


■ 


INTERNAL  LOAD  (KN) 
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Figure  E.3  Ram  Friction  Calibration  Curve 


APPENDIX  F 


Problems  encountered  with  the  laboratory  eau i omen t 
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The  high  pressures  used  in  the  laboratory  program  required 
consideration  to  be  given  to  all  sources  of  error  resulting 
from  leaks  in  the  apparatus,  membranes  and  lines.  Poulos 
(1964)  discusses  various  sources  of  leaks  and  these  are 
summarized  below. 

F low  of  Fluid  through  Membranes 

Flow  through  membranes  can  be  caused  by  hydraulic 
pressure  gradients  ,  vapour  pressure  and  mole  fraction 
gradients.  The  rate  of  flow  of  water  is  directly 
proportional  to  the  hydraulic  pressure  difference.  Soaking 
of  membranes  in  silicone  oil  prior  to  the  test  resulted  in 
equal  or  higher  permeabilities  than  those  soaked  in  water. 

This  problem  of  flow  through  membranes  is  of  more 
concern  in  long  term  tests  than  short  term  tests.  The 
duration  of  tests  in  this  program  is  less  than  5  hours  and 
this  source  of  leaks  is  of  little  significance. 

Leakage  past  Bi ndi ngs 

This  source  of  error  is  again  of  more  importance  for 
long  term  tests.  Smooth  pedestals  with  a  single  0-ring  on 
each  side  of  the  porous  disc  resulted  in  leakage  of 
approximately  0.1  ml/day,  which  is  negligible  for  short  term 
tests . 

The  membrane  clamping  system  was  modified  in  this 
program  to  reduce  the  number  of  bindings  and  0-rings 
required.  The  author  was  concerned  more  with  large  leaks 


. 
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developing  (due  to  puncturing  of  membranes  or  poor  seals) 
than  these  small  leaks  observed  by  Poulos. 

Leakage  from  Va 1 ves  and  F i 1 1 i nqs 

Poulos  showed  that  Klinger  valves  are  generally 
unsuitable  for  any  type  of  laboratory  work.  The  Klinger 
valves  that  were  in  use  with  the  pressure  systems  used  by 
this  author  were  observed  to  leak  noticeably  at  relatively 
low  pressures  and  were  thus  replaced  by  circle  seal  valves. 
These  valves  showed  no  visible  signs  of  leakage.  The 
majority  of  valves  used  in  the  plumbing  systems  are  Swagelok 
no-volume  change  valves  and  these  have  also  performed 
satisfactory . 

Leaks  through  valves  and  fittings  will  be  critical  in 
long  term  tests  as  small  leaks  will  not  be  noticeable  due  to 
rapid  evaporation  of  the  fluid.  However,  careful 
calibration  of  the  volume  change  apparatus  can  accomodate 
all  these  types  of  leaks.  The  apparatus  in  use  in  this 
program  showed  no  measureable  volume  changes  occur ing  over  a 
24  hour  period. 

•  Ruptur i nq  of  Membranes 

This  was  the  biggest  problem  faced  by  the  author  during 
his  laboratory  program.  At  pressures  exceeding  5.6  MPa,  the 
internal  membrane  tends  to  pinch  between  the  top  of  the 
sample  and  the  load  cap  or  porous  stone.  A  thin  Teflon  ring 
placed  over  such  joints  proved  partially  successful,  but  it 
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was  found  that  the  membrane  developed  small  holes  in  the 
vicinity  of  this  ring.  The  teflon  ring  was  replaced  by  a 
Neoprene  rubber  ring  which  has  performed  satisfactory  thus 
far . 

The  internal  membrane  has  been  the  dominating  factor  in 
the  many  tests  that  have  had  to  be  abandoned  due  to 
development  of  large  leaks.  The  thickness  of  the  membrane 
has  been  increased  but  this  has  not  solved  the  problem 
completely.  Any  imperfections  in  the  membrane  will  result 
in  a  hole  developing  due  to  sand  grain  penetration  at  high 
pressure.  A  redesign  of  this  membrane  (and  probably  of  the 
load  cap)  will  be  required  if  a  testing  program  is  to  be 
conducted  successfully  on  a  regular  basis. 

The  external  membrane  has  performed  satisfactory.  Only 
one  hole  developed  and  this  was  probably  due  to  a  weakness 
along  the  seam. 

Pressure  Systems  used  i n  the  Laboratory 

The  hydrogen  pressure  system  used  for  the  internal 
pressure  operated  well  and  maintained  pressure  to  within  a 
few  kPa.  The  biggest  drawback  of  this  system  is  the  size  of 
its  accumulator  and  reservoir  (1100  ml).  The  total  volume 
changes  occurring  during  the  entire  test  are  large,  and  most 
of  this  volume  is  lost  from  the  reservoir  if  a  leak  occurs 
in  the  membrane  as  the  pressure  system  has  to  be  shut  down. 
The  reservoir  has  no  indicator  for  the  oil/water  interface 
and  it  is  thus  possible  to  pump  oil  through  the  system. 
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It  is  suggested  that  a  larger  accumul ator /reserved  r 
system  be  provided  for  future  testing  and  that  some  means  of 
monitoring  the  oil/water  interface  be  installed. 

The  other  pressure  system  initially  used  (  Wykeham) 
proved  unsatisfactory  as  it  did  not  maintain  pressure  very 
well.  Replacement  parts  were  not  available  as  it  is  now 
outdated.  This  system  was  finally  replaced  by  a  constant 
pressure  pump  which  has  an  accumulator  built  into  the 
system.  This  pressure  system  has  also  operated  exceedingly 
well  and  has  a  reservoir  which  can  be  filled  manually  during 
the  test. 
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Grain  size  distr ibut ions  and  compress ibi 1 i tv  data 
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'Figure  G.1  Grain  size  distribution  for  the  Oedometer  sample 
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GRAIN  SIZE  (MM) 


Figure  G.2  Grain  size  distribution  for  Test  #A-1 
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Figure  G.3  Grain  size  distribution  for  Test  #A-2 
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Figure  G.4  Grain  size  distribution  for  Test  #B-1 
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Figure  G.5  Grain  size  distribution  for  Test  #B-2 
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Figure  G.6  Grain  size  distribution  for  Test  #C-1 
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Figure  G.7  Grain  size  distribution  for  Test  #C-2 
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Figure  G.8  Compressibility  test,  Sample  #A-1 


Figure  G.9  Compressibility  test,  Sample  #A-2 
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Figure  G.10  Compressibi 1 i ty  test,  Sample  #B-1 
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Figure  G.11  Compressibility  test,  Sample  #B-2 
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Figure  G.12  Compressibility  test,  Sample  #C-1 
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Figure  G.13  Compressibility  test,  Sample  #C-2 
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